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Effects  of  Mixtures  on  Liquid  and  Solid  Fragment  Size  Distributions 


Culbert  B.  Laney^ 

Engility  Corp.,  8211  Terminal  Rd,  Lorton,  VA  22079  U.S.A. 


Abstract:  This  paper  discusses  the  effects  of  mixing  diverse  fragment  populations.  Such  mixing 
naturally  occurs  when  the  tensile  strain  rates  that  drive  fragmentation  vary  with  time  and/or 
space.  Understanding  such  mixtures  requires  understanding  the  ingredients,  i.e.,  simple 
fragmentation  driven  by  steady  uniform  strain  rates  in  a  fixed  dimension.  Theoretical  proofs  are 
given  showing  that  simple  fragmentation  -  as  well  as  its  inverse,  simple  coagulation  -  leads  to 
an  infinite  but  discrete  set  of  universal  size  distributions.  This  conclusion  is  a  middle  ground 
between  previous  studies  which  found  one  universal,  two  universals,  or  an  infinite  continuum  of 
non-universals.  This  paper  finds  that  complex  multi-stage  fragmentation  events  may  often  be 
treated,  in  aggregate,  as  virtual  simple  single-stage  fragmentation  events. 

Key  Words:  Fragmentation,  coagulation,  universal  size  distributions,  Rosin-Rammler 
distributions,  Weibull  distributions.  Gamma  distributions,  Mott-Linfoot  distributions,  Grady- 
Kipp  distributions,  Simmons  root  normal  distribution,  maximum  entropy  theory,  Grady’s  law, 
mining,  milling,  orbital  debris,  hypervelocity  impact,  weapon  effects,  atomization  and  sprays 

1.  Introduction 

Fragmentation  occurs  in  many  different  applications  including:  rock  blasting  for  mining  and 
excavation;  materials  processing  such  as  milling,  grinding,  and  crushing;  orbital  debris  from 
artificial  satellites  including  collisional  cascades;  hypervelocity  impacts  and  collisional  evolution 
of  asteroids  and  meteorites;  cratering,  ejecta,  dust,  and  debris  from  high-explosive  and  nuclear 
weapon  detonations;  cased  munitions  including  sympathetic  detonations;  debris  from  high-speed 
projectile  penetration;  bubble-bearing  magmatic  ejecta  from  volcanic  eruptions;  burn  rates  in 
fluidized-bed  combustors;  and  liquid  atomization  and  sprays. 

Researchers  have  found  that  manifestly-different  fragmentation  events  may  obtain  the  same  or 
similar  outcomes.  For  example,  Domokos  et.  al.  (2015)  report  that  “fragmenting  various  types  of 
materials  by  slowly  proceeding  weathering  and  by  rapid  breakup  due  to  explosion  and 
hammering”  produces  fragments  size  distributions  and  shapes  with  “astonishing  universality  . . . 
in  spite  of  the  strongly  different  cracking  mechanisms  governing  fragment  formation  in  dynamic 
fragmentation  initiated  by  explosion  and  hammering  and  in  weathering  induced  spallation.”  For 
another  example,  in  a  literature  survey.  Fancy  (2016)  found  a  number  of  cases  where  liquids  and 
solids  obtained  identical  or  nearly-identical  fragment  size  distributions.  This  was  true  despite  the 
fact  that  liquids  and  solids  fragment  in  entirely  different  ways.  In  particular,  liquids  typically 
fragment  via  stable  and  unstable  flows  followed  by  secondary  breakup  and  coagulation,  while 
solids  typically  fragment  via  crack  propagation. 

Other  researchers  have  found  that  seemingly-similar  fragmentation  events  may  obtain  radically 
different  outcomes.  For  example,  in  a  comprehensive  survey  of  “size  distributions  of  asteroid 
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families  . . .  [believed  to  be]  remnants  of  parent  bodies  that  broke  into  fragments,”  Parker  et.  al. 
(2008)  found  size  distributions  with  power  law  exponents  ranging  continuously  from  —0.1  to 
—1.04.  On  this  basis,  they  concluded  that  the  “size  distribution  varies  significantly  [and]  ... 
cannot  be  described  by  a  universal  function.”  In  addition,  they  found  that  many  “size 
distributions  . . .  display  a  well-defined  change  of  slope  and  can  [only]  be  modeled  as  a  ‘broken’ 
double  power  law.” 

This  paper  attempts  to  explain  these  divergent  findings  by  distinguishing  between  simple  events 
that  produce  a  single  type  of  fragment  and  compound  events  that  produce  mixtures  of  different 
fragments  formed  under  different  conditions  at  different  times,  in  different  places,  and  in 
different  dimensions.  Following  Laney  (2016),  this  paper  draws  freely  on  results  previously 
obtained  for  various  solids,  liquids,  and  hypothetical  materials,  such  as  those  studied  in  Discrete 
Element  Methods  (DEMs),  Molecular  Dynamics  (MD),  and  other  types  of  models.  Thus  this 
paper  benefits  from  the  progress  made  in  a  wide  variety  of  different  communities. 

2.  Size  Distributions 

Suppose  fragment  sizes  can  be  described  as  follows: 

X-c,D'  (1) 

where  i  and  c,  are  constants.  Notice  that  €■=  1  for  /  =  1  and  /  >  3 .  Also  notice  that  c_-=  1/C; .  The 
three  most  common  choices  of  X  are  as  follows: 


i  =  1 

X 

=  D 

diameter 

i  =  m  —  l 

X 

surface  area 

i  =  m 

X 

=  M  =  pV’ 

mass 

where  1  <  m  <  3  is  the  fragment  dimension,  p  is  density,  and  "V  is  volume.  Eess  common  but 
still  important  choices  of  X  are  as  follows: 

i  =  -\  X  =  \l D  inverse  diameter 

i  =  -(m  - 1)  X  =1/^  inverse  surface  area 

i  =  -m  X  =1/ M  inverse  mass 

In  general,  liquids  produce  spherical  or  nearly-spherical  fragments;  thus  liquid  fragments  are 
usually  characterized  hy  X  =  D  (or  D  ^ ).  By  contrast,  in  general,  solids  produce  rough  irregular 
fragments.  Thus  solid  fragments  are  usually  characterized  by  X  =  M  (or  M  ' ). 


2 


Because  it  is  difficult  measure  directly,  the  surface  area  ^  is  rarely  used  to  characterize  either 
liquid  or  solid  fragments.  However,  as  discussed  in  Appendix  A,  <51^  is  a  natural  choice,  because 
the  work  required  to  create  fragments  is  proportional  to 

For  example,  for  spherical  fragments  such  as  liquid  droplets: 

m  =  3;  Cj^=pnl6 

For  another  example,  for  solids  such  as  rock,  Domokos  et.  al.  (2015)  found  that  “above  a 
characteristic  size  the  overall  shape  of  all  fragments  can  be  well  approximated  by  rectangles  of 
dimensions  1:  1.56  :  2.32.”  Let  D  be  the  middle  dimension  and  let: 

AR^  =1/1.56  =  0.641  and  AR^  =2.32/1.56  =  1.487 

be  the  aspect  ratios.  Then: 

m  =  3;  C2  =  2(A/?j +A/?2 +A/?2ARj)  =  6.163;  03=  AR2ARj/?  =  0.953/? 

For  comparison,  perfectly  square  fragments  obtain: 

m  =  3  ;  C2  =  6. ;  c^=  p 

As  a  final  example,  Wittel  et.  al.  (2006)  measured  “regular  isotopic”  solid  shell  fragments  to  find 
m  =  2  and  “needle-like”  plate  glass  fragments  to  find  m  =  1.5,  which  “implies  that  fragments 
have  a  self-affine  character,  meaning  that  the  larger  they  are,  the  more  elongated  they  get.” 

Because  researchers  rarely  specify  m,  it  must  usually  be  estimated  from  experimental 
photographs  or  written  descriptions  of  fragment  shapes.  More  specifically,  if  the  aspect  ratio 
does  not  depend  on  fragment  size,  and  if  the  fragments  retain  two  parallel  surfaces  from  the 
original  body,  then  m  —  2.  Alternatively,  if  the  aspect  ratio  does  not  depend  on  fragment  size, 
and  if  the  fragments  do  not  retain  any  surfaces  from  the  original  body,  then  m  =  3  . 

Fragment  size  distributions  may  be  described  by  f(X),  i.e.,  the  number  fraction  of  fragments  in  a 
range  dX  centered  on  X  divided  by  dX.  In  standard  probability  theory,  f(X)  is  called  the 
probability  density  function  (PDF).  Notice  that  f(X)  is  always  non-negative  such  that: 


\f(X)dX=l  (2a) 

0 

Alternatively,  fragment  size  distributions  may  be  described  by  /j.(X),i.e.,  the  Y  fraction  of 
fragments  in  a  range  dX  centered  on  X  divided  by  dX  where  Y  =  .  Notice  that  fy  (2f )  is 

always  non-negative  such  that: 
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(2b) 


\fy{X)dX=\ 

0 


Also  notice  that: 


fy{X) 


Yf(X) 

00 

jYf(X)dX 

0 


(3) 


3.  Average  Fragment  Sizes 

Consider  the  following  averages  based  on  diameter,  surface  area,  mass,  and  so  forth: 


X,.,=\xf(X)dX  (4) 

0 


\XYf(X)dX 

Xy..,  =  J  Xf,(X)dX  =  iL -  (5) 

°  \Yf(X)dX 

0 


where  X  =c.D‘ 


and  Y  -  CjD^ .  The  best  known  examples  are  as  follows: 


D,.,=\mD)dD  (6) 

0 

00 

DM..,=lDf„{D)dD  (7) 

0 

which  are  known  as  the  count  mean  diameter  (CMD)  and  the  mass  mean  diameter  (MMD), 
respectively.  Less  common  but  still  important,  consider  the  following  averages  based  on  inverse 
diameter,  inverse  surface  area,  inverse  mass,  and  so  forth: 
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1 

^=1 

avg  0 


fjX) 

X 


dX 


(8) 
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i  X 


’jr/ix) 


dX 


]Yf{X)dX 


The  best  known  example  is  as  follows: 


1 


/«(£>) 


D  •’  D 

‘^U  avg  0 


dD 


or  equivalently: 


D[ 


M  avg  CO 


/m(^) 

D 


dD 


which  is  called  the  Sauter  mass  mean  diameter  (SMD). 
Consider  the  following  ratios  of  average  fragment  diameters: 


Qx  = 


D 


X  avg 


D 


avg 


avg 


Dl 


avg 


R  = 


D' 


avg 


R.  = 


D 


X  avg 


D' 


X  avg 


(9) 


(10a) 


(10b) 


(11) 


(12) 


(13) 


For  Equation  (13),  the  best  known  example  is  as  seen  in,  e.g.,  Simmons  (1977),  Wu  et.  al. 
(1991),  Chou  &  Faeth  (1998),  and  Sallam  et.  al.  (2006). 

Finally,  consider  the  following  ratio  of  average  fragment  sizes: 


ref 


CyD 


(14) 


where  may  be  X^^^,  Xy^^^ 


,or  x: 


Yavg  and  may  be  Dy^^^,  or  Dy 


Y  avg  • 
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4.  Rosin-Rammler  Size  Distributions 


4.1  Introduction 


Rosin-Rammler  size  distributions  may  be  defined  as  follows: 


/(£>)  = 


AD 


ref 


exp 


^  D  ^ 


y^refj 


(15a) 


where  a.  A,  I,  and  n  are  parameters  and  is  a  reference  diameter.  If  Y  =  CjD^ ,  Equation  (15a) 
can  be  rewritten  as  follows: 


fy(D)- 
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^Y^ref 
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exp 


t  A" 

D 


(15b) 


If  X  =c.D‘ ,  Equations  (15a)  and  (15b)  can  be  rewritten  as  follows: 
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y  ^ref  J 


exp 


-\a 


y  ^ref  J 


fyiX)  =  -. 


iAyXref 
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exp 
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(16a) 


(16b) 


By  Equation  (2): 


l+\ 

fz+n 

FT  H  "  ^ 

r 

y  n  J 

1 

/+7+1 

r/+j+q 

— 

a  «  r 

y  n  ) 

(17) 


(18) 


Eaney  (2015b)  suggested  three  self- similarity  conditions: 

1.  When  a  size  distribution  is  expressed  in  terms  of  a  given  average  size,  it  should  actually 
obtain  that  average  size. 

2.  When  a  size  distribution  is  expressed  in  terms  of  a  given  average  size,  and  that  average 
size  changes,  key  free  parameters  in  the  size  distribution  should  remain  the  same. 

3.  When  a  size  distribution  is  expressed  in  terms  of  a  given  average  size,  and  that  average 
size  changes,  the  size  distribution  should  stay  the  same. 
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Table  1  shows  parameters  for  Rosin-Rammler  size  distributions.  This  table  assume  that  the  first 
self-similarity  condition  applies  and  that  Z  =  c^D^ . 


Table  1.  Parameters  for  Rosin-Rammler  size  distributions  assuming  the  first  self-similarity  condition  applies. 
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The  second  self- similarity  condition  is  true  for  Rosin-Rammler  size  distributions  if: 


l  =  k=l2  (19a) 

n  =  n,  =^2  (19b) 


where  the  subscripts  1  and  2  refer  to  any  two  reference  diameters.  As  it  turns  out,  Equation  (19) 
also  ensures  the  third  self-similarity  condition,  at  least  for  the  averages  described  in  Section  3. 
To  prove  this,  notice  that  the  second  and  third  self-similarity  conditions  are  both  true  if: 
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This  is  true  if: 


“  ^2^ref2 

Multiply  the  first  equation  by  the  second  equation  taken  to  the  {l  +  \)l  n  power  to  obtain: 


/+1  /+1 
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By  Equation  (17),  this  is  true  because: 


l+l  l+l 
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Table  2  shows  expressions  for  the  ratios  of  fragment  diameters  assuming  the  self-similarity 
conditions  apply.  As  before,  Table  2  assumes  X  -  . 


Table  2.  Ratios  of  average  fragment  diameters  for  Rosin-Rammler  size  distributions  assuming  the  self- 
similarity  conditions  apply. _ 
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The  following  subsections  focus  on  three  specific  classes  of  Rosin-Rammler  size  distributions 
that  have  proven  physical  meanings. 

4.2  Weibull  Size  Distributions 

As  defined  here,  a  Weibull  size  distribution  can  be  written  in  the  following  form: 

const.  ejq3 [- const,  ] 

where  x  is  some  measure  of  fragment  size.  As  one  special  case,  consider  Rosin-Rammler  size 
distributions  with  I  -  n-l  and  a  >  0 : 


or  equivalently: 
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^  D  ^ 


n-l 


AD 


ref 


y^refj 


exp 


/  A" 

D 
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(20a) 


(20b) 


where,  as  usual,  Y  -  CjD^ .  Equations  (20a)  and  (20b)  can  be  rewritten  as  follows: 
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Table  3.  Parameters  for  Type  II  Weibull  size  distributions  assuming  the  first  self-similarity  condition  applies. 


Table  4  shows  expressions  ratios  of  fragment  diameters  assuming  the  self-similarity  condition 
applies.  As  before,  X  =  . 


Table  4.  Ratios  of  average  fragment  diameters  for  Type  II  Weibull  size  distributions  assuming  the  self¬ 
similarity  conditions  apply. 


9 


As  another  special  case,  consider  Rosin-Rammler  size  distributions  with  I  =  n  —  j  — I  and  a  >  0 : 
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exp 
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K^ref  J 
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or  equivalently: 
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exp 
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(24b) 


where,  as  usual,  Y  =  .  Equations  (24a)  and  (24b)  can  be  rewritten  as  follows: 
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where,  as  usual,  X  =  c-D‘  and  Y  -  CjD^ .  By  Equations  (17)  and  (18): 
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(27) 


Table  5  shows  results  for  other  parameters  assuming  that  /  =  n-  j  -I  and  that  the  first  self¬ 
similarity  condition  applies.  As  before,  Z  -c^D'‘  . 
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Table  5.  Parameters  for  Type  I  Weibull  size  distributions  assuming  the  first  self-similarity  condition  applies. 


With  the  parameter  choices  shown  in  Tables  3  and  5,  Weibull  size  distributions  satisfy  all  three 
self-similarity  conditions.  In  particular,  the  key  parameter  n  does  not  depend  on  . 

Not  all  combinations  of  i,  j,  k,  and  n  are  allowed.  For  example,  suppose  and 

j  =  k  =  0 .  Then  Tables  3  and  5  give  the  following: 


a  =  r 


f  1  y 

i+- 
V  nj 


which  is  valid  if  and  only  if  n<-l  or  n>0.  For  another  example,  suppose  and 

j  =  k  =  0 .  Then  Tables  3  and  5  give  the  following: 


r 

a  =  Y  1 

V 


n) 


which  is  valid  if  and  only  if  n  <  0  or  n  >  1 . 


Notice  that,  for  n<0,  Weibull  size  distribution  become  power  law  size  distributions  when 
D  »  .  This  treatment  assumes  that  power  law  size  distributions  are,  in  fact,  the  large  size 

limit  of  Weibull  size  distributions  with  n<0;  see  also  Laney  (2015b). 

As  seen  in  Laney  (2015b),  the  most  common  Weibull  forms  can  be  designated  as  Types  I  and  II. 
More  specifically.  Type  I  Weibull  size  distributions  can  be  written  as  follows: 
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/m(^)  = 
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ref 
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and  Type  II  Weibull  size  distributions  can  be  written  as  follows: 
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The  Type  I  is  obtained  from  Equations  (24)  and  (25)  with  X  —Y  —  M  or,  equivalently, 

i  =  j  =  m.  The  Type  II  is  obtained  from  Equations  (20)  and  (21)  with  X  -  M  ox,  equivalently, 

i-m. 


Weibull  size  distributions  with  n  >  0  are  often  associated  with  ductile  materials  while  those  with 
n  <  0  are  often  associated  with  brittle  shatter-prone  materials.  However,  as  seen  later,  any 
material  -  including  simple  liquids  like  water  -  can  fragment  with  either  n  >  0  or  n<0.  Instead 
of  ductile  vs.  brittle,  which  can  change  with  strain  rate,  this  treatment  distinguishes  n>0  from 
n<0  based  on  fragment  geometry. 

As  seen  in  Table  6,  there  is  an  association  between  Type  II  Weibull  size  distributions  and 
geometry,  a.k.a.,  perfect  packing  theory.  More  specifically.  Type  II  Weibull  distributions  with 
n  >  0  are  associated  with  cases  where  adjacent  fragments  meet  along  extended  lines  or  surfaces, 
so  that  the  parts  to  fit  together  like  jigsaw  puzzle  pieces,  e.g.,  tessellations,  subdivisions. 
Similarly,  Type  II  Weibull  size  distributions  with  n  <  0  are  associated  with  cases  where  adjacent 
fragments  meet  only  at  a  single  point,  e.g.,  fractals,  Apollonian  sphere  packings. 
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Table  6.  Examples  of  Type  II  Weibull  size  distributions  derived  from  geometric  proofs.  Cowan  (2010)  and 
Gilvarry  (1961)  are  exact  while  the  others  are  numerical  approximations.  For  negative  n,  the  original  results 
are  power  laws. _ _ _ _ 


nim 

n 

m 

Description 

Reference 

1 

2 

2 

Various  random  subdivisions  by 
straight  lines  based  on  Poisson 
distributions. 

Grady  &  Kipp  (1985) 
Cowan  (2010) 

-1/2  to-1 

-1  to  -2 

2 

Various  random  2D  circle,  triangle, 
square,  and  diamond  packings.  The 
most  common  value  of  n  is - 1.3. 

Aste  (1996) 

-2/3 

-2 

3 

Random  subdivision  by  arbitrary 
surfaces,  lines,  and  points  based  on 
mutually-independent  Poisson 
distributions  assuming  the  fragments 
are  “geometrically  similar.” 

Gilvarry  (1961) 
Gilvarry  & 
Bergstrom  (1961) 

-2/3  to  -  1 

-2  to  -  3 

3 

Various  random  3D  Apollonian 
sphere  packings.  The  most  common 
value  of  n  is - 2.5 

Eind  et.  al.  (2008) 

As  shown  by  Laney  (2015b),  Type  I,  Type  II,  and  other  types  of  Weibull  size  distribution  mimic 
each  other,  except  for  very  large  or  very  small  fragments.  As  a  result,  without  exact  solutions,  it 
is  hard  to  distinguish  the  real  size  distributions  from  their  close  approximations.  Fortunately, 
there  is  at  least  one  exact  solution.  More  specifically,  for  certain  kinds  of  linear  subdivisions, 
Cowan  (2010)  proved  that  a  Type  II  Weibull  distribution  -  and  not  any  other  type  of  Weibull 
size  distribution  -  is  an  exact  geometric  solution. 

More  importantly.  Brown  &  Woheltz  (1995)  gave  an  exact  semi-geometric  proof  for  Type  II 
Weibull  distributions  in  general.  First,  they  introduced  a  simplified  form  of  the  Population 
Balance  Equation  (PBE)  involving  a  “single-event  particle  distribution  function;”  see  also,  e.g., 
Hamilton  et  al.  (2003)  and  Kostoglou  &  Karabelas  (2004).  Next,  based  on  geometric  theory,  they 
argued  that  the  single-event  particle  distribution  function  is  a  power  law;  see  also,  e.g.,  Turcotte 
(1986).  Einally,  they  proved  that  Type  II  Weibull  distributions  are  the  exact  analytic  solutions  of 
their  simplified  Population  Balance  Equation. 

Based  on  the  available  evidence,  it  will  be  assumed  that  Type  II  Weibull  size  distributions 
uniquely  represent  fragment  geometries. 

4.3  Gamma  Size  Distributions 

As  defined  here,  a  Gamma  size  distribution  can  be  written  in  the  following  form: 

const. e?q3[-|Z7|x”J 

where  x  is  some  measure  of  fragment  size.  As  one  special  case,  consider  Rosin-Rammler  size 
distributions  with  I  -  a  - 1 : 
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Equations  (30a)  and  (30b)  can  be  rewritten  as  follows: 
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ref 
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where,  as  usual,  X  -  c^D'  and  Y  =  .  These  will  be  called  Type  II  Gamma  size  distributions; 

see  also  Laney  (2015b).  Notice  that  earlier  treatments,  including  Laney  (2015b),  were  limited  to 
positive  a.  However,  this  treatment  allows  both  negative  and  positive  a. 

Section  4. 1  gives  expressions  for  n ,  A ,  Aj,  and  as  functions  of  a  and  .  For  example, 

see  Table  7.  As  shown  by  Laney  (2015b),  with  the  parameter  choices  given  in  Table  7,  Type  II 
Gamma  size  distributions  satisfy  the  first  and  third  self-similarity  conditions,  but  not  the  second 
self-similarity  condition.  In  other  words,  the  key  parameter  a  depends  on  the  choice  of  . 
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Table  7  gives  implicit  expressions  for  n  as  a  function  of  a.  In  general,  these  must  be  solved 
iteratively.  However,  if  a  >  0  andD^^^  -  ,  there  is  a  simple  analytical  solution: 


n  —  1 


which  follows  from  r(a  +  1)  =  ar(a) .  Similarly,  if  a  <  0  and  ,  there  is  a  simple 

analytical  solution: 

n  —  —l 


As  another  special  case,  consider  Rosin-Rammler  size  distributions  with  I  =  iaS^  -1 : 
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or  equivalently: 
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exp 
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If  X  —M  ,  these  will  be  called  Type  IV  Gamma  size  distributions.  This  is  a  short  but  opaque 
definition.  Laney  (2015b)  gives  a  longer  but  clearer  definition  of  Type  IV  Gamma  size 
distributions. 

More  generally,  every  X  is  associated  with  a  different  type  of  Gamma  size  distribution. 
However,  regardless  of  type.  Gamma  size  distributions  satisfy  the  first  and  third  self-similarity 
conditions,  but  not  the  second  self-similarity  condition.  In  other  words,  the  key  parameter  a 
always  depends  on  . 

Similar  to  Weibull  size  distributions,  there  is  a  strong  association  between  Gamma  size 
distributions  and  geometrical  theory.  In  fact,  the  research  literature  suggests  the  three  possible 
geometrical  interpretations: 

1.  Kiang  (1966)  argued  that  Type  IV  Gamma  size  distributions  approximate  fragment 
size  distributions.  As  evidence,  he  showed  that  Type  IV  Gamma  size  distributions 
approximate  the  results  of  Voronoi  tessellation.  The  association  between  Gamma  size 
distributions  and  Voronoi  tessellation  has  been  confirmed  by  numerous  subsequent 
researchers,  e.g.,  Ferenc  &  Neda  (2007).  However,  the  association  between  Voronoi 
tessellation  and  fragmentation  has  not  been  confirmed.  To  the  contrary,  Grady  &  Kipp 
(1985)  observed  that  while  “Kiang  has  very  strongly  stated  his  belief  that  the  Voronoi 
method  is  at  the  foundation  of  random  fragmentation,  ...  we  personally  feel  that  the 
Voronoi  construction  introduces  too  much  correlation  among  neighboring  points  to 


15 


provide  a  good  description  of  fragmentation.”  Having  said  this,  there  is  limited  evidence 
that  Voronoi  tessellation  represents  cases  where  both  fragmentation  and  its  inverse, 
coagulation,  occur;  see,  e.g..  Appendix  A.l  in  Laney  (2016). 

2.  Melzak  (1953)  argued  that  Type  IV  Gamma  size  distributions  are  the  natural  initial 
conditions  for  coagulation.  More  specifically,  he  showed  that  the  Smoluchowski 
coagulation  equation  has  an  analytical  solution  if  the  initial  size  distribution  is  a  Type  IV 
Gamma  size  distribution;  see  also  Friedlander  &  Wang  (1966),  Scott  (1968),  and 
Lindblad  (2005,  2007). 

3.  Villermaux  et.  al.  (2004)  and  Marmottant  &  Villermaux  (2004)  argued  that  Type  II 
Gamma  size  distributions  are  the  natural  final  conditions  for  the  Smoluchowski 
coagulation  equation.  The  proof  is  similar  to  that  given  by  Brown  &  Woheltz  (1995),  as 
discussed  above. 

This  treatment  will  adopt  the  third  and  most  modern  interpretation.  In  fact,  the  proof  given  by 
Villermaux  et.  al.  (2004)  and  Marmottant  &  Villermaux  (2004)  appears  be  equally  valid  for 
Gamma  size  distributions  associated  with  diameter  (Type  II),  inverse  diameter,  surface  area, 
inverse  surface  area,  mass  (Type  IV),  inverse  mass,  and  so  forth.  In  other  words,  based  on  the 
available  evidence,  it  will  be  assumed  that  various  types  of  Gamma  size  distributions  represent 
coagulation  geometries. 

Recent  research  has  shown  that,  in  some  cases,  fragmentation  is  inextricably  linked  with 
coagulation.  For  example,  Villermaux  et.  al.  (2004)  argued  that  a  common  liquid  “fragmentation 
mechanism  ...,  somewhat  surprisingly,  consists  of  a  coalescence  process.”  Similarly,  Marmottant 
&  Villermaux  (2004)  note  that  liquid  “blobs,  just  before  breakup,  interact,  and  that  the 
interaction  is  of  a  coalescence  or  aggregation  type.”  Less  commonly,  solids  can  experience 
something  like  coagulation  due  to  surface  deformation  such  as  occurs  in  explosive  consolidation 
of  powders,  surface  melting,  surface  wetness,  electrostatic  forces,  magnetic  forces,  chemical 
forces,  gravitational  forces,  van  der  Waals  forces,  and  so  forth;  see,  e.g.,  Brilliantov  (2015), 
Vledouts  et.  al.  (2015,  2016a). 

4.4  Exponential  Power  Law  Size  Distributions 

As  defined  here,  an  exponential  power  law  size  distribution  can  be  written  in  the  following  form: 

const,  exp [- const,  J 

where  x  is  some  measure  of  fragment  size.  As  one  special  case,  consider  Rosin-Rammler  size 
distributions  with  I  —  i-1  and  a>0: 


f(X)  = 


iAX 


-exp 


ref 


-a 


V  ) 


(33a) 
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c  i  ^  Y 

f  =  ^ 


exp  -  a 


(33b) 


where,  as  usual,  X  =  c.D‘  and  Y  =  CjD^ .  Section  4.1  gives  expressions  for  a.  A,  Ay  and 
as  functions  of  n  and  X  . .  For  example,  by  Equations  (17)  and  (18): 


1  --  i 

A  =  —  a  ” r  - 
\n\  V  n 


A^=Xa  -  r( 


As  another  special  case,  consider  Rosin-Rammler  size  distributions  with  I  =  i  —  j  —  \  and  a  >  0 : 


/F(^)=./y  exp  -a 

ref  [  X  ref 


(36a) 


or  equivalently: 


/(X)  =  ^^  ^ 

iAXref[X,,, 


exp  -  a 


(36b) 


where,  as  usual,  X  =  c^D‘  and  Y  =  c-D^ .  As  before.  Section  4.1  gives  expressions  for  a  ,  A , 
Ay  and  as  functions  of  n  and  X  . .  For  example,  by  Equations  (17)  and  (18): 


1  i  -  i 

A  =  ^a  ”  r 

\n\  V  n 


A  =  A  a''  t(  - 


With  the  parameter  choices  shown  in  Tables  1  and  2,  exponential  power  law  size  distributions 
satisfy  all  three  self-similarity  conditions.  In  other  words,  the  key  parameter  n  does  not  depend 
on 
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The  following  is  a  generalization  of  a  maximum  entropy  proof  given  by  Cousin  et.  al.  (1996). 
Entropy  is  defined  as  follows: 


CO 

5  =  -const,  j  fiX)\n[fiX)]dX  (39) 

0 

Based  on  Equation  (2),  two  possible  constraints  are  as  follows: 

CO 

\fiX)dX=l  (40) 

0 


\fy(X)dX=l  (41) 

0 


where,  as  usual,  X  -  c.D'  and  Y  -CjD^ .  Using  Lagrange  multipliers  to  maximize  entropy 
subject  to  these  two  constraints: 


/(X)  =  const. exp(-const.X^^‘)  (42) 

Notice  that  Equation  (42)  is  the  same  as  Equation  (33a)  if  n  =  j.  Also  notice  that  Equation  (42) 
does  not  automatically  ensure  Equations  (40)  and  (41).  Rather,  the  two  constants  must  be  chosen 
specifically  to  ensure  Equations  (40)  and  (41).  While  Equation  (40)  is  usually  enforced. 

Equation  (41)  is  rarely  enforced. 

As  one  example,  suppose  that: 


X=Y  =  M 


Then  Equations  (40)  and  (41)  become: 


00 

^f(M)dM=l 

0 

00 

\Mf{M)dM=M^^^ 

0 


In  addition.  Equation  (42)  becomes: 

f{M)  =  const. exp(-constAf) 


(43) 


(44) 


(45a) 


or  equivalently: 
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f{D)  =consti)'"  ^  exp  (-const  Z)'") 


(45b) 


Notice  that  Grady  &  Kipp  (1985)  obtained  Equation  (45a)  for  m  —  2  while  Li  &  Tan  kin  (1987) 
obtained  Equation  (45b)  for  m-3  . 

As  alternatives  to  Equations  (40)  and  (41),  two  possible  constraints  are  as  follows: 

CO 

\f(X)dX=l  (46) 

0 


J/i  (X)dX=l 

0  Y 

Using  Lagrange  multipliers  to  maximize  entropy  subject  to  these  two  constraints: 

f(X)  =  const. exp(-const.X^^^') 

Notice  that  this  is  the  same  as  Equation  (33a)  if  n  =  —j.  As  one  example,  suppose  that: 

X  =Y  =  M 

Then  Equations  (46)  and  (47)  become: 

00 

j/(M)ZM  =1 

0 


J 


/(M) 

M 


dM 


1 


M', 


avg 


In  addition.  Equation  (48)  becomes: 


r 

/(M)  =  const  .exp 

V 


const.  ^ 

M  y 


Einally,  suppose  the  two  constraints  are  as  follows: 


\f,{X)dX=\ 

0 


(47) 


(48) 


(49) 


(50) 


(51) 


(52) 
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(53) 


\f^^AX)dX=\ 

0 

where,  as  usual,  X  =  c,D' ,  Y  -  CjD^ ,  and  Z  =  .  Using  Lagrange  multipliers  to  maximize 

entropy  subject  to  these  two  constraints: 

/j,(X)  =  const. exp(-const.X“^^')  (54) 

This  is  the  same  as  Equation  (36a)  if  n  =  ±k. 

Suppose  that  more  than  two  constraints  are  imposed.  For  example,  Ahmadi  &  Sellens  (1993) 
imposed  four  constraints  as  follows: 


oo 

\fiD)dD  =  l 

0 

CO  *^1 

I  fi,D(D)dD  =  const. J  —  /(Z))JZ)  =  1 

0  Q  ^ 


J  (D)dD  =  const  .J  D^f(D)dD  =  1 

0  0 


J  /^3  (D)dD  =  const.J  D^f(D)dD  =  1 

0  0 


(55) 

(56) 

(57) 

(58) 


and  obtained: 


f 

/(D)  =  const  .exp 

V 


const. 

D 


const.D^  -  const  D^ 

y 


(59) 


In  this  case,  and  in  general,  multiple  constraints  lead  to  exponentials  of  linear  combinations  of 
powers  of  D.  Notice  that  Equation  (59)  does  not  automatically  ensure  Equations  (55)  to  (58). 
Rather,  the  four  constants  must  be  chosen  specifically  to  ensure  Equations  (55)  to  (58),  which 
requires  solving  a  4x4  system  of  equations. 

For  more  information  on  maximum  entropy  theory  see,  e.g.,  the  surveys  by  Englman  (1996), 
Babinsky  &  Sojka  (2002),  Dumouchel  (2009),  and  Dechelette  et.  al.  (2011).  Based  on  the 
available  evidence,  it  will  be  assumed  that  exponential  power  laws  satisfy  maximum  entropy 
constraints  and  are  equally  relevant  to  fragmentation  and  coagulation. 
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5.  Universal  Size  Distributions 


As  defined  here,  universal  size  distributions  for  fragmentation  are  those  size  distributions  which 
are  both  Type  II  Weibull  and  exponential  power  laws.  As  seen  in  Section  4.2,  Type  II  Weibull 
size  distributions  satisfy  constraints  imposed  by  fragmentation  geometry  while  exponential 
power  law  size  distributions  satisfy  constraints  imposed  by  maximum  entropy  theory.  Thus 
universal  size  distributions  for  fragmentation  satisfy  both  geometric  and  maximum  entropy 
constraints  simultaneously. 

In  addition  to  exponential  power  laws,  maximum  entropy  theory  yields  all  sorts  of  other  size 
distributions,  e.g..  Equation  (59).  However,  only  exponential  power  laws  intersect  Type  II 
Weibull  size  distributions.  For  earlier  discussions  of  the  intersections  between  Type  II  Weibull 
and  exponential  power  law  size  distributions  -  or,  equivalently,  the  intersections  between 
fragment  geometry  and  maximum  entropy  theory  -  see,  for  example,  Grady  &  Kipp  (1985), 
Brown  &  Wohletz  (1995),  and  Kumar  &  Kumaran  (2005). 

By  Equations  (20)  and  (33),  Type  II  Weibull  size  distributions  are  the  same  as  exponential  power 
law  size  distributions  if  and  only  if  /  =  ( - 1  =  n  - 1  or: 


n  -  i  (60) 

Thus  the  problem  becomes  one  of  selecting  X  =c^D\  seen  in  Section  2,  the  most  common 
and  physical  choices  of  X  are  as  follows: 

i  =  ±\  or  X  -  D-' 

i  =  ±(m  —  1)  or  X  = 

i  -  ±m  or  X  = 

Type  II  Weibull  size  distributions  with  n-i  -  ±l,  ±  (m  - 1) ,  and  ±  m  will  be  called  the 
fundamental  size  distributions. 

Table  8  summarizes  the  fundamental  size  distributions  for  integer  m.  For  m  =  2 ,  D  =  ^  and 
thus  there  are  only  two  fundamental  size  distributions.  Similarly,  for  m  =  1,  <%  =  const. and 
D  =  M  and  thus  there  is  only  one  fundamental  size  distribution. 


Table  8.  Fundamental  size  distributions  for  integer  m. 


m-3 

m-2 

m  =  1 

n  -±1 

n!  m  —  ±1/3 

Mott  &  Einfoot  (1943) 

n! m-  ±1/2 

Mott  &  Einfoot  (1943) 

n!  m  —  ±\ 

Eineau  (1936) 

n  — ±2 

n!  m  —  ±2/3 

Gilvarry  (1961) 

n!  m  —  ±\ 

Grady  &  Kipp  (1985) 

n  -  ±3 

n!  m  —  ±\ 

Ei  &  Tan  kin  (1987) 
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Notice  that  all  of  the  references  given  in  Table  8  concern  positive  fundamentals,  except  for 
Gilvarry  (1961),  which  concerns  a  negative  fundamental.  In  addition,  all  of  the  references  in 
Table  8  provide  full  theoretical  proofs,  except  for  Mott  &  Linfoot  (1943),  who  provide  only  a 
partial  proof.  These  earlier  proofs  generally  differ  from  the  one  given  here. 

Some  rows  and  columns  in  Table  8  have  been  known  for  decades.  For  example.  Tucker  et.  al. 
(1965)  described  the  first  row  of  Table  8.  More  specifically,  they  noted  that  a  Type  II  Weibull 
size  distribution  with  n! m  —  \  “applies  for  random  or  semi-random  breakup  of  a  bar,” n! m  —  \l  2 
“applies  for  breakup  of  a  thin  metal  shell  and  is  observed  to  hold  for  two-dimensional  breakup,” 
and  n!  m  —  H  3  “applies  for  a  very  thick  metal  shell  and  is  observed  to  hold  for  three- 
dimensional  breakup  [and]  . . .  may  be  expected  to  hold  for  the  smallest  fragments.”  For  another 
example,  Gilvarry  (1961)  described  something  similar  to  the  first  column  of  Table  8.  More 
specifically,  for  m  =  3  ,  he  suggested  Type  I  Weibull  size  distributions  with  n  =  1,  2,  and  3. 
However,  he  is  best  remembered  for  deriving  a  Type  II  Weibull  size  distribution  with  n  =  -2. 

Turning  from  fragmentation  to  coagulation,  universal  size  distributions  for  coagulation  are  those 
size  distributions  which  are  both  Gamma  and  exponential  power  laws.  As  seen  in  Section  4.3, 
Gamma  size  distributions  satisfy  constraints  imposed  by  coagulation  geometry  while  exponential 
power  law  size  distributions  satisfy  constraints  imposed  by  maximum  entropy  theory.  Thus 
universal  size  distributions  for  coagulation  satisfy  both  geometric  and  maximum  entropy 
constraints  simultaneously. 

By  Equations  (31)  and  (33),  Type  II  Gamma  size  distributions  are  the  same  as  exponential  power 
law  size  distributions  if  and  only  if  I  =  ( - 1  =  a  - 1  or: 


a  -  i 


As  for  fragmentation,  the  three  most  physical  choices  of  i  are  1,  m-1 ,  and  m  .  However,  using 
Table  7,  it  can  be  shown  that  there  is  a  solution  for  n  if  and  only  if: 


a  —  i  —  n  —  \ 


and: 


(62) 


(63) 


Using  the  same  approach  on  other  types  of  Gamma  size  distributions,  the  universal  size 
distributions  for  coagulation  are  the  same  as  those  for  fragmentation.  The  only  difference  is  that 
coagulation,  unlike  fragmentation,  restricts  the  reference  size;  see  Table  9. 

Table  9.  Examples  of  restrictions  on  reference  sizes  implied  by  universal  size  distributions  for  coagulation. 


n-\ 

D 

ref  avg 

n--\ 

D 

ref  avg 

n  —  m-1 

n  —  -(m-1) 

'^ref  avg 

n  =  m 

n  =  -m 

M  =M' 

ref  avg 
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The  restrictions  given  in  Table  9  resolve  an  inherent  ambiguity  with  Gamma  size  distributions. 
More  specifically,  Gamma  size  distributions  no  longer  violate  the  second  self-similarity 
condition  because  the  key  parameter  a  no  longer  depends  on  an  arbitrary  reference  size  but  rather 
on  a  single  clearly-defined  reference  size. 

The  remainder  of  this  paper  is  limited  to  integer  m  .  If  the  fundamentals  are  n  =  ±1, ...,  ±  m ,  then 
the  overtones  are  Type  II  Weibull  size  distributions  where: 

n  -  ±(m  + 1),  ±  (m  -I-  2),  ±  (m  -l-  3), . . . 

In  addition,  the  semitones  or  halftones  are  Type  II  Weibull  size  distributions  where: 

n  =  ±l/2,±3/2,  ±5/2,... 


The  following  sections  argue  that,  for  m  —  3,  the  universal  size  distributions  consist  of  the 
fundamentals,  overtones,  and  semitones.  In  other  words: 

n  =  ±I/2,±I,  ±3/2,  ±2,  ±5/2,  ±3,  ±7/2, ... 


or  equivalently: 


n  1112  5,7 

m  6  3  2  3  6  6 


Because  these  values  are  integer  multiples  of  one  sixth,  this  may  be  called  the  law  of  sixths.  The 
law  of  sixths  also  seems  to  apply  for  m  =  1  and  2.  In  other  words,  only  the  classification  of  the 
universals  changes  with  m,  e.g.,  nl m-112  is  a  fundamental  for  m  =  2  but  a  semitone  for 
m  —  3. 


6.  Simple  Fragmentation 

As  defined  here,  simple  fragmentation  occurs  when  the  tensile  strain  rate  s  and  fragment 
dimension  m  are  constants.  This  assumes  that  the  tensile  strain  rate  drives  fragmentation  and  that 
it  can  be  expressed  as  a  scalar  rather  a  full  tensor.  As  seen  in  Section  A. 9,  examples  include 
radially-expanding  thin  cylindrical  shells,  radially-expanding  thin  spherical  shells,  radially- 
expanding  thin  rings,  stretched  straight  rods,  stretched  linear  filaments,  and  so  forth. 

Appendix  A  shows  that  simple  fragmentation  leads  to  simple  algebraic  expressions  for  average 
fragments  sizes.  In  addition,  it  will  be  argued,  simple  fragmentation  leads  to  universal  size 
distributions  about  those  average  fragment  sizes.  Mott  &  Linfoot  (1943)  were  the  first  to  make 
the  association  between  universal  size  distributions  and  simple  algebraic  expressions  for  average 
fragment  sizes. 

As  noted  in  the  introduction,  many  communities  study  fragmentation.  However,  only  two 
communities  routinely  study  simple  fragmentation,  namely,  the  conventional  weapons  effects 
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community  in  terms  of  cased  munitions,  and  the  liquid  atomization  and  sprays  community  in 
terms  of  round  jets,  planar  jets,  air  bubbles,  rings,  ligaments,  and  so  forth. 


Table  10  gives  examples  of  size  distributions  taken  from  the  research  literature  on  cased 
munitions  where  n! m,  n  ,  and  m  refer  to  parameters  in  Type  II  Weibull  size  distributions.  In 
these  examples,  the  cylindrical  shells  obtain  fundamentals  while  the  rings  obtain  overtones.  As 
described  in  Appendix  A,  simple  fragmentation  only  occurs  for  sufficiently  thin  rings  and  shells. 
This  not  very  common;  see  sections  7.3.2  and  8.5  for  examples  involving  thick  shells. 


Table  10a.  Fragment  size  distributions  for  tests  of  solid  uniformly-expanding  thin  rings  with  m=l. 


nl  m 

n 

m 

Reference 

Notes 

2 

2 

I 

Zhou  et.  al.  (2006) 

Various  tests  found  in  a  literature  survey 

4 

4 

I 

Moxnes  &  Bprve  (2015) 

25mm-thick  steel  rings  cut  from  warheads 

Table  10b.  Fragment  size  distributions  for  tests  of  solid  uniformly-expanding  thin  cylindrical  shells  with  m=2. 


nl  m 

n 

m 

Reference 

Notes 

1/2 

1 

2 

Bess  (1975) 

10  to  17mm-thick  SAE  1050  steel 

It 

It 

It 

Arnold  &  Rottenkolber  (2008) 

2mm-thick  steel;  n/m  =  0.48, 0.55 

-1/2 

-1 

2 

Bannikova  et.  al.  (2014) 

2mm-thick  alumina;  n/m- -0.47 ± 0.02 

Table  10c.  Fragment  size  distributions  for  tests  of  solid  uniformly-expanding  thin  cylindrical  shells  with  m=3. 


n/m 

n 

m 

Reference 

Notes 

1/3 

1 

3 

Arnold  &  Rottenkolber  (2008) 

6mm-thick  steel;  n/m  =  0.29, 0.31, 0.35 

2/3 

2 

3 

Grady  et.  al.  (2001) 

8%-thick  heat-treated  steel;  n/m  =  0.67 

Table  1 1  gives  examples  of  size  distributions  taken  from  the  research  literature  on  atomization 
and  sprays  where  n/m,  n  ,  and  m  refer  to  the  parameters  in  Type  II  Weibull  size  distributions 
as  converted  from  the  Type  I  root  normal  and  Type  II  Gamma  size  distributions  originally  given 
in  most  of  the  references;  see  Laney  (2015a,  2016).  Notice  that  all  the  tests  obtained 
fundamentals  except  the  last,  which  obtained  an  overtone. 


Table  11.  Fragment  size  distributions  for  various  tests  involving  liquids  with  m=3. 


n/m 

n 

m 

Reference 

Notes  (%  of  tests  obtaining  distribution) 

1/3 

1 

3 

Spielbauer  et.  al.  (1989) 

Sheet  jets  (60%) 

2/3 

2 

3 

Wuet.  al.  (1991) 

Cylindrical  jets  (16.6%) 

2/3 

2 

3 

Marmottant  &  Villermaux  (2004) 

Stretched  ligaments  (33.3%) 

1 

3 

3 

Bremond  &  Villermaux  (2006) 

Sheet  jets  (16.6%) 

1 

3 

3 

Lhuissier  &  Villermaux  (2012) 

Spherical  shells  (“air  bubbles”)  (100%) 

-1 

-3 

3 

Lhuissier  et.  al.  (2013) 

Radial  ligaments  from  drop  impact  (100%) 

2 

6 

3 

Chou  &  Faeth  (1998) 

Semi-spherical  bags  (50%) 

In  summary,  the  available  experimental  evidence  indicates  that  simple  fragmentation  always 
obtains  either  a  fundamental  or  an  overtone. 
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7.  Compound  Simple  Fragmentation 

As  defined  here,  compound  simple  fragmentation  is  any  non-simple  case  that  results  in  a  single 
overall  Type  II  Weibull  size  distribution.  Two  common  outcomes  are  universal  size  distributions, 
as  discussed  earlier,  and  compound  size  distributions.  For  universal  size  distributions,  n! m  tends 
to  remain  constant  regardless  of  the  exact  conditions.  When  it  does  change,  it  jumps  from  one 
discrete  value  to  another.  By  contrast,  for  compound  size  distributions,  n! m  changes 
continuously  as  the  exact  conditions  change.  Sections  7.1  to  7.4  give  examples  of  universal  and 
compound  size  distributions  based  on  a  literature  survey.  Section  7.5  gives  a  basic  theory, 
involving  mixing  and  the  integral  mean  value  theorem,  for  why  universal  and  compound  size 
distributions  occur  in  compound  simple  fragmentation. 

7.1  One  Dimension  {m  =  1) 

One-dimensional  objects  such  as  bars,  rings,  rods,  etc.  can  be  loaded  in  various  different  ways 
including:  axially  in  tension;  axially  in  compression;  and  laterally.  As  seen  Section  6,  simple 
axial  tension  results  in  simple  fragmentation,  i.e.,  fundamentals  and  overtones  with 
n! m  —  ±1,  ±  2,  ±3, . . . 

As  an  example  of  axial  compression.  Gladden  et.  al.  (2005)  studied  “dynamic  buckling  and 
fragmentation  of  slender  rods  axially  impacted  by  a  projectile  ...  for  a  range  of  materials 
including  teflon,  dry  pasta,  glass,  and  steel.”  They  obtained  an  expression  for  average  fragment 
size  similar  to  those  given  in  Appendix  A.  However,  while  the  results  are  similar,  the  proofs  are 
different.  More  specifically.  Appendix  A  assumes  fragmentation  is  governed  by  the  tensile 
component  of  the  usual  conservation  laws,  while  Gladden  et.  al.  (2005)  assumes  fragmentation  is 
driven  by  unstable  oscillatory  bending  modes.  Gladden  et.  al.  (2005)  give  experimental  size 
distributions  that  appear  to  be,  at  best,  piecewise  compound  simple  and  possibly  completely  non¬ 
simple. 

As  an  example  of  lateral  loading,  Ching  et.  al.  (2000)  studied  impact  fragmentation  of  long  thin 
glass  rods  with  lengths  of  2m  and  diameters  of  2mm.  When  the  glass  rods  were  dropped 
horizontally  onto  the  ground  from  heights  of  1.2  or  2.0m,  they  found  that:  (a.)  the  fragments 
formed  in  one-dimension  and  (b.)  a  single  power  law  (equivalently,  a  single  Type  II  Weibull  size 
distribution  with  n<0)  provided  a  good  fit  for  entire  fragment  size  distribution.  When  the  glass 
rods  were  dropped  from  larger  heights  up  to  23.6m,  they  found  that  the  smallest  fragments 
formed  in  three  dimensions  and  obeyed  one  power  law,  while  the  largest  fragments  formed  in 
one  dimension  and  obeyed  another  power  law.  Of  the  seven  power  law  indices  reported,  five 
approximately  agree  with  the  law  of  sixths.  More  specifically,  two  have  nl m  ^  —113  and  three 
have  n/mx— 1/2.  This  indicates  that  when  fragments  form  in  a  full  three-dimensional  space, 
they  tend  to  obtain  fragment  size  distributions  typical  of  three-dimensions,  regardless  of  the  fact 
that  m  —  1. 
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7.2  Two  Dimensions  {m  =  2) 


7.2.1  Case  1  {nim  =  -1/3) 

To  simulate  space  debris,  Wittel  et.  al.  (2004,  2005,  2006)  studied  fragmentation  due  to  small 
hydrogen  gas  explosions  inside  eggshells.  In  all  cases  tested,  they  obtained  a  power  law 
(equivalently  a  Type  II  Weibull  size  distribution)  with  m-2  and  n/m^  -1/ 3 .  According  to 
Figure  I  in  Wittel  et.  al.  (2005),  “the  ignition  took  place  at  the  top  of  the  egg,”  causing  the  top  of 
the  eggshell  to  fragment  first  and  relieving  overpressure  on  the  bottom  of  the  eggshell.  Since  the 
top  experienced  a  higher  strain  rate  than  the  bottom,  as  opposed  to  the  uniform  strain  rates 
characteristic  of  simple  fragmentation,  this  is  an  example  of  compound  simple  fragmentation. 

As  another  way  to  simulate  space  debris,  Wittel  et.  al.  (2004,  2005,  2006)  studied  “intact 
eggshells  . . .  catapulted  onto  the  ground  at  a  high  speed  using  a  simple  setup  of  rubber  bands.”  In 
all  cases  tested,  they  obtained  m-2  and  n/mx—\/3.  According  to  Figure  2  in  Wittel  et.  al. 
(2005),  when  an  “egg  [shell]  hits  the  ground  in  the  direction  of  its  longer  axis  . . .  [the]  egg  [shell] 
suffers  gradual  collapse  as  it  moves  forward,  making  the  impact  process  relatively  long  ...”  The 
timing  between  sequential  experimental  photographs  indicates  that  the  eggshell  front  experiences 
higher  strain  rates  than  the  back.  In  addition,  the  eggshell  experiences  primarily  lateral  loading 
on  the  front  and  back,  and  compressive  loading  on  the  sides,  as  opposed  to  tensile  loading  in 
explosive  cases.  The  fact  that  such  radically  different  loading  conditions  -  both  explosive  and 
impact  -  always  obtained  n/m  «  -1/3  implies  that  this  is  a  true  universal  size  distribution. 

In  addition  to  tests,  Wittel  et.  al.  (2004,  2005,  2006)  modeled  egg  shell  fragmentation  using  a 
Discrete  Element  Method  (DEM)  for  a  hypothetical  material  composed  of  “pointlike  material 
elements  [with]  . .  .bonds  between  nodes  . . .  assumed  to  be  springs  having  linear  elastic  behavior 
up  to  failure.”  This  model  obtained  n/m-  -0.35  ±  0.03  for  impact  fragmentation  and 
film-  -0.55  ±  0.03  for  explosive  fragmentation,  which  approximately  follows  the  law  of  sixths. 

7.2.2  Case  2  {nIm  =  -5/6) 

In  an  early  study  of  space  debris,  Bess  (1975)  used  a  light  gas  gun  to  fire  a  1.65g  steel  cylinder  at 
3.0km/s  and  a  0.37g  aluminum  cylinder  at  4.5  km/s  into  a  simulated  spacecraft  wall  consisting  of 
an  “insulated  fiberglass  wall  in  back  of  which  were  a  number  of  electronic  boxes  containing 
resistors,  capacitors,  etc.”  Despite  the  complex  geometry  and  the  variety  of  different  wall 
materials ,  the  two  tests  each  obtained  a  single  fragment  size  distribution.  More  specifically,  the 
wall  fragments  obeyed  a  power  law  with  nlm~  -0.80  for  the  steel  test  and  n/m  «  -0.84 for  the 
aluminum  test.  Bess  described  the  wall  fragments  as  “irregularly  shaped,  fiat  plates,”  which 
implies  m-2.  Based  on  these  and  four  subsequent  tests,  the  NASA  Standard  Breakup  Model 
assumes  n  =  -1.71  for  collisions  involving  satellites  or,  equivalently,  n/m  «  -0.855 ;  see,  e.g., 
Johnson  et.  al.  (2001),  Sakuraba  et.  al.  (2008).  This  approximately  obeys  the  law  of  sixths. 

7.2.3  Case  3  (n/m  =  -5/4) 

The  previous  cases  approximately  obeyed  the  law  of  sixths,  which  assumes  the  same  set  of 
universal  size  distributions  occur  regardless  of  whether  m  =  1 , 2,  or  3.  However,  this  is  only  true 
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when  fragmentation  occurs  in  a  full  three-dimensional  space.  As  an  example  that  does  not  obey 
the  law  of  sixths,  Kun  &  Herrmann  (1999)  modeled  collisions  between  two-dimensional  disks 
using  a  Molecular  Dynamics  (MD)  model.  This  model  treats  each  disk  as  a  hypothetical  material 
composed  of  “unbreakable,  undeformable  ...  randomly  shaped  convex  polygons”  connected  by 
breakable  elastic  beams.  Because  this  model  confined  fragmentation  to  two  dimensions,  it 
obtained  a  two-  rather  than  a  three-dimensional  semitone.  More  specifically,  all  but  the  largest 
fragments  obtained  a  Type  II  Weibull  size  distribution  with  n/m  —  —1.27  ±0.05  «  — 5/4  or, 
equivalently,  nx—5/2.  This  is  a  semitone  for  m-2,  i.e.,  it  is  halfway  between  the  fundamental 
n  —  -2  and  the  overtone  n  —  -3. 

7.3  Three  Dimensions  (m  =  3) 

This  section  considers  six  cases  that  obey  the  law  of  sixths,  as  summarized  in  Table  12,  and  one 
case  that  violates  the  law  of  sixths. 


Table  12.  Universal  size  distributions  for  various  compound  simple  fragmentation  tests  with  m=3. 


nl  m 

m 

Reference 

Notes 

1 

-0.20 

+0.06 

3 

Timar  et.  al.  (2010) 

Impact  fragmentation  of  small  (4mm) 
polypropylene  spheres. 

-0.15 

-0.17 

-0.20 

3 

Oddershedde  et.  al.  (1993) 

Impact  fragmentation  of  thin  (4mm) 
gypsum  plates. 

2 

-1/3 

3 

Edwards  &  Deal  (2011) 

Explosive  fragmentation  of  9.25mm-thick 
aluminum  alloy  cylindrical  shell. 

3 

-0.48 

3 

Oddershedde  et.  al.  (1993) 

Impact  fragmentation  of  frozen  7.5cm 
cubes  of  gypsum,  stearic  paraffin,  soap, 
and  potato 

-0.49 

+0.18 

3 

Onose  &  Fujiwara  (2004) 

Porous  gypsum  spheres  impacted  by  nylon 
spheres  fired  from  a  light  gas  gun. 

4 

-2/3 

3 

Gilvarry  &  Bergstrom  (1961) 
Yamakoshi  et.  al.  (1981) 

Inaoka  &  Takayasu  (1996,  99) 
Hooper  (2012) 

Impact  fragmentation  of  glass  spheres, 
asteroids,  meteorites,  planetesimals,  and 
consolidated  aluminum  powder  spheres. 

5 

-0.83 

+0.03 

3 

Seebaugh  (1977) 

Schoutens  (1979) 

Ivanov  et.  al.  (1983) 

O’Keefe  &  Ahrens  (1985,  87) 
Turcotte  (1986) 

Rubbleized  and  ejected  rock  from  high- 
seed  impact,  high-explosive  and  nuclear 
tests  such  as  Elat  Top,  Pile  Driver,  and 
Danny  Boy. 

0.84^ 

3 

Vledouts  et.  al.  (2016b) 

Explosive  fragmentation  of  thin  roughly- 
spherical  water- glycerol  shells 

6 

-7/6 

3 

Kaminski  &  Jaupert  (1998) 

Explosive  volcanic  eruptions  of  bubble¬ 
bearing  magma. 

-1.15 

+0.08 

3 

Durand  &  Soulard  (2012,  13) 

Shock-loaded  metal  melts  resulting  in 
planar  jets. 

1.  Converted  from  Type  II  Gamma  size  distributions  as  in  Laney  (2015b,  2016) 
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7.3.1  Case  1  {nim  =  -1/6) 


This  subsection  concerns  Case  1  in  Table  12.  As  one  example,  to  model  “industrial  processes 
related  to  mining  and  ore  processing,”  Timar  et.  al.  (2010)  studied  isotactic  polypropylene 
spheres  with  diameters  of  4mm  impacting  a  hard  wall.  Based  on  experimental  results  for  impact 
velocities  of  40,  60,  and  75  m/s,  they  suggested  a  Type  II  Weibull  size  distribution  (equivalently 
a  power  law)  with  ri/ m  -  -0.2  +  0.06 .  As  another  example,  Oddershedde  et.  al.  (1993)  studied 
4mm-thick  gypsum  plates  “thrown  horizontally  onto  the  floor  [or]  . . .  thrown  vertically  down” 
and  obtained  n/m  =  -0.15  to  -0.20.  In  all  these  cases, n/m  = -1/6  to  within  the  experimental 
error.  These  examples  concern  impact  fragmentation.  Here  and  in  general,  impact  fragmentation 
tends  to  produce  highly  unsteady  and  non-uniform  tensile  strain  rates.  A  possible  exception  is 
when  the  flat  face  of  a  rectangular  object  impacts  a  flat  surface. 

7.3.2  Case  2  {nIm  =  -1/3) 

This  subsection  concerns  Case  2  in  Table  12.  Edwards  &  Deal  (2011)  studied  explosive-filled 
7075-T6  cylindrical  aluminum  alloy  shells.  The  as-received  (“AR”)  shells  were  heat-treated  in 
two  different  ways  (“OA”  and  “ST”)  to  reduce  their  strength  and  hardness.  The  outer  radius  of 
the  shells  varied  while  all  other  parameters  were  kept  fixed.  More  specifically,  all  shells  had  a 
length  of  33m,  an  inner  diameter  of  6.5mm,  and  were  fragmented  by  detonating  a  35mm-long 
1.05g  granulated  pellet  of  tetryl.  Figure  1  shows  a  typical  result. 


M(g) 

Figure  1.  Type  II  Weibull  distribution  vs.  test  data  for  explosive-filled  aluminum  alloy  cylindrical  shells  with 
an  outer  radius  of  25mm  as  reported  by  Edwards  &  Deal  (2011). 

Table  13  summarizes  the  best-fit  Type  II  Weibull  size  distributions  (equivalently  power  laws). 
Notice  that  all  size  distributions  obey  the  law  of  sixths.  However,  the  values  given  in  the  shaded 
cells  are  somewhat  uncertain  due  to  the  coarse  binning  used  in  the  tests.  More  specifically. 
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Edwards  &  Deal  (201 1)  used  only  12  size  bins  with  fixed  limits,  meaning  that  as  few  as  four  size 
bins  were  actually  occupied  in  any  given  test. 


Table  13.  Fragment  size  distributions  based  on  curve  fits  to  test  data  reported  by  Edwards  &  Deal  (2011). 


Case 

Strengt 

h  (MPa) 

Outer  Diameter  (mm)  [Explosive-to-Metal 

1  Mass  Ratio] 

Yield 

Tensile 

33  [0.0130] 

25  [0.0231] 

20  [0.0382] 

15  [0.0752] 

“OA” 

274+41 

343+1.4 

n/m  =  -2/3 

n/m  =  -1/3 

n/m  =  -1/2 

n/m  =  —1/2 

“ST” 

401+15 

551+22 

n/m  =  -2/3 

n/m  =  -1/3 

n/m  =  -1/2 

n/m  =  —1/2 

“AR” 

519+13 

578+12 

n/m  =  -2/3 

n/m  =  -1/2 

n/m  =  -1/2 

n/m  =  —1/2 

In  addition  to  the  cases  shown  in  Table  13,  Edwards  &  Deal  (201 1)  also  tested  thicker  shells  with 
outer  radii  of  42mm,  50mm,  and  100mm.  In  the  42mm  case,  the  strongest  shell  produced  eight 
fragments,  the  middling  shell  produced  three  fragments,  and  weakest  shell  cracked  but  did  not 
fragment.  The  50mm  and  100mm  shells  did  not  fragment  regardless  of  material  strength. 

Due  to  the  thick  shell  walls,  the  explosive  charge-to-shell  mass  ratio  in  Edwards  &  Deal  (201 1) 
was  much  smaller  than  in  most  similar  tests.  Eor  example,  the  ratio  was  2  to  1 1  times  smaller 
than  in  Mock  &  Holt  (1981,  1983).  This  may  explain  why  Mock  &  Holt  (1981,  1983)  obtained 
n/m  >  0,  as  discussed  in  Section  8.5,  while  Edwards  &  Deal  obtained  n/m<0. 

7.3.3  Case  3  {n/m  =  -1/2) 

This  subsection  concerns  Case  3  in  Table  12.  To  simulate  high-speed  collisions  involving 
asteroids  and  planetesimals,  Onose  &  Eujiwara  (2004)  studied  64%-porous  gypsum  spheres 
impacted  by  nylon  spheres  fired  from  a  two-stage  light  gas  gun.  They  found 
film  —  —0.49  ± 0. 1 8  w  1  / 2  for  early-time  fragments  “ejected  conically  within  a  few  msec  after 
the  impact”  and  n!  m  —  —\  .49  ±  0.09  «  3  /  2  for  late-time  fragments  “consisting  of  hundreds  of 
slow,  small  fragments  ejected  almost  perpendicular  to  the  target.”  Separating  fragments  in  time 
is  rare  -  most  experiments  mix  early-  and  late-time  fragments  together.  The  early-time  results 
agree  with  those  obtained  by  Oddershedde  et.  al.  (1993)  when  “stearic  paraffin,  soap,  and  potato 
were  frozen  in  liquid  nitrogen  before  being  fragmented  by  a  hammer.”  Oddershedde  et.  al. 

(1993)  concluded  that  “within  the  accuracy  of  the  determination  of  the  exponents  no  significant 
dependence  on  material  was  observed.”  This  insensitivity  implies  that  n! m  —  -1/2  is  a  true 
universal  size  distribution. 

7.3.4  Case  4  {nim  =  -2/3) 

This  subsection  concerns  Case  4  in  Table  12.  As  an  early  example,  Gilvarry  &  Bergstrom  (1961, 

1 962)  obtained  n  /  m  =  — 2  /  3  for  fragmentation  of  small  glass  spheres  impacted  by  a  piston.  They 
included  gelatin  “to  prevent  the  fragments  from  colliding  with  the  walls  of  the  retaining 
chamber;  otherwise,  they  will  undergo  secondary  fracture  by  virtue  of  their  large  kinetic  energy.” 
In  this  regard,  they  noted  that  “. . .  experimental  [results]  ...  for  single  fracture  were  virtually 
nonexistent  in  the  literature”  at  that  time. 

As  another  example,  Yamakoshi  et.  al.  (1981)  obtained  n/m-  —2/ 3  for  asteroid  break  up, 
planetesimal  fragmentation,  and  fragmentation  of  iron  meteorites;  see  also  Grady  &  Kipp  (1987). 
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As  a  third  example,  based  on  a  literature  survey  and  a  three-dimensional  impact  fracture  model, 
Inaoka  &  Takayasu  (1996,  1999)  concluded  that  n/m  =  -2/3  “seems  to  be  a  very  universal 
value.”  As  a  fourth  example.  Hooper  (2012)  obtained  n! m  —  —2l for  “fragmentation  of 
brittle,  granular  aluminum  spheres  following  high  velocity  impact  (0. 5-2.0  km/s)  on  thin  steel 
plates.”  More  specifically,  he  studied  “spherical  porous  aluminum  projectiles  2.51  cm  in 
diameter  . . .  cut  from  a  cylinder  of  Valimet  H-2  aluminum  powder  which  had  been  isostatically 
pressed  into  a  monolith.”  As  a  final  example,  Domokos  et.  al.  (2015)  found  n! m  —  —'ll 3  in 
seven  different  fragmentation  experiments  involving  rocks  and  minerals  such  as  limestone, 
dolomite,  and  gypsum.  The  fragmentation  was  done  in  substantially  different  ways  including 
gradual  weathering,  hammering,  and  explosions.  Domokos  et.  al.  (2015)  observe  that  “it  is 
astonishing  that  the  same  universality  class  is  recovered”  in  all  cases. 

o 

Astrom  et.  al.  (2004)  noted  that  “within  the  mining  engineering  community,  . . .  [this]  scale- 
invariant  fragment-size  distribution  has  long  been  known.”  However,  modern  mainstream 
mining  engineering  practice  typically  involves  an  infinite  continuum  of  compound  size 
distributions  rather  than  universal  size  distributions;  see,  e.g.,  Cunningham  (2005). 

7.3.5  Case  5  (n/m  =  ±5/6) 

This  subsection  concerns  Case  5  in  Table  12.  Seebaugh  (1977),  Schoutens  (1979),  Ivanov  et.  al. 
(1983),  O’Keefe  &Ahrens  (1985,  1987),  and  Turcotte  (1986)  surveyed  test  results  for  rock 
fragmentation  due  to  high-velocity  impacts,  large  chemical  explosions,  and  nuclear  tests.  In  most 
cases,  they  found  power  laws  with  a  ~  0.5  where: 

a  =  n  +  m. 

In  other  words,  in  most  cases  they  found  Type  II  Weibull  size  distributions  with  n/m  «  -5/6  .  In 
one  case,  O’Keefe  &  Ahrens  (1985)  found  that  the  smallest  fragments  obtained  nimx  -0.633 
while  the  largest  fragments  obtained  nimx  -0.827.  Notice  that  the  former  value  is 
approximately  a  fundamental  while  the  latter  value  is  approximately  a  semitone. 

For  an  essentially  unbounded  fragmenting  object  such  as  the  earth,  the  majority  of  fragments  will 
be  produced  in  peripheral  regions  where  the  fragmentation  forces  barely  exceed  the  material 
strength,  while  only  a  minority  of  fragments  will  be  produced  in  the  core  regions  where  the 
fragmentation  forces  greatly  exceed  the  material  strength.  This  is  true  regardless  of  how  large  the 
fragmenting  forces  may  be.  This  means  that  negative  size  distributions  characteristic  of  brittle 
fragmentation  will  tend  to  be  observed  rather  than  the  positive  distributions  characteristic  of 
ductile  fragmentation,  even  in  the  most  extreme  cases,  such  as  high-yield  nuclear  weapon 
detonations. 

O’Keefe  &  Ahrens  (1987)  speculated  that  “additional  comminution”  including  “grinding  and 
crushing  . . .  occurs  after  the  initial  shock  has  passed”  leading  to  a  “broader  size  distribution  in 
impact  ejecta  as  compared  to  that  obtained  in  simple  brittle  failure  experiments.”  While  trying  to 
verify  a  similar  hypothesis,  Kaminski  &  Jaupart  (1998)  repeatedly  dropped  a  piston  onto  samples 
of  rhyolitic  pumice.  Regardless  of  the  number  of  piston  strikes,  they  found  that  the  fragment 
“population  always  keeps  the  same  power  law  exponent  of  n  =  -2.6+0. 1.”  Assuming  m  =  3,  this 
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approximately  corresponds  to  n/m  —  — 5/6 .  Finally,  out  of  frustration,  they  put  the  fragments 
through  a  rotating  grinder,  either  with  or  without  a  small  steel  ball.  Not  surprisingly,  this  reduced 
the  fragment  size  spread  and  increased  the  absolute  value  of  the  power  law  exponent.  The  fact 
that  n/m  =  -5/6  appears  in  so  many  different  kinds  of  tests,  and  persists  through  multiple 
stages  of  fragmentation,  implies  that  it  is  a  true  universal  value. 

As  a  final  example,  Vledouts  et.  al.  (2016b)  used  rapid  exothermic  gaseous  reactions  to  drive 
expansion  of  thin,  roughly-spherical  liquid  shells.  They  found  that  hydrodynamic  instability 
leads  to  “holes  [that]  grow  in  radius  by  capillary  retraction,  accumulating  the  liquid  of  the  shell 
into  their  rim,  which  merge  with  the  rim  of  the  neighboring  hole,  finally  leaving  the  shell  as  a 
web  of  ligaments,  which  break  by  capillary  instability;”  see  also  their  Figure  4.  The  various 
circular  rims  and  ligaments  probably  each  experiences  simple  fragmentation.  However,  because 
the  different  rims  and  ligaments  form  at  different  times  and  places,  and  thus  experience  different 
strain  rates,  the  net  result  is  compound  simple  fragmentation  with  n/m  «  5/6.  While 
n/m  =  5/6  may  seem  to  be  very  different  than  n/m  =  -5/6 ,  the  results  given  in  Sections  8.2 
and  8.3  imply  that  they  are  actually  closely  related. 

7.3.6  Case  6  (n/m  =  -7/6) 

This  subsection  concerns  Case  6  in  Table  12.  Kaminski  &  Jaupert  (1998)  note  that  “explosive 
volcanic  eruptions  eject  a  large  mass  of  magma  fragments,  ranging  from  meter-sized  blocks  to 
micron-sized  ash  particles  dispersed  at  high  altitudes  in  the  atmosphere.”  In  a  literature  survey, 
Kaminski  &  Jaupert  (1998)  found  eleven  instances  of  powers  laws  with  n  «  -3  and  six  instances 
of  power  laws  with  n  «  -3.5 .  Assuming  m-3 ,  these  correspond  to  n/m  «  -1  and 
n/m  «  -7/6,  respectively,  where  the  former  is  a  fundamental  and  the  latter  is  a  semitone. 

Although  bubble -bearing  liquid  magma  might  seem  to  be  classically  ductile,  Kaminski  &  Jaupert 
(1998)  found  only  power  laws,  which  are  typically  indicative  of  brittle  fragmentation.  Indeed, 
Clarke  (2013)  notes  that  “high  magma  acceleration  rates  can  generate  large  strain  rates  that  may 
result  in  a  rheological  transition  from  ductile  to  brittle  behavior  . . .  causing  the  magma  to  behave 
in  a  brittle  fashion,  and  leading  to  fragmentation.” 

As  another  example,  Durand  &  Soulard  (2012,  2013)  studied  fragmentation  of  shock  loaded 
metal  melts.  More  specifically,  the  specified  loading  conditions  created  three  sheet  jets: 

“During  the  expansion  of  the  sheets,  holes/cradcs,  randomly  distributed,  begin  to 
appear...  Once  the  void  has  percolated,  a  complementary  network  of  ligaments  of  liquid 
metal  between  the  pores  appears  ...In  the  late  times,  the  oblique  views  of  the  whole 
systems  show  that  the  ligaments  previously  created  in  the  three  sheets  have  broken  up 
and  have  reached  a  final  spherical  form.  ” 

The  various  liquid  metal  ligaments  probably  each  experience  simple  fragmentation.  However, 
because  they  fragment  at  different  strain  rates,  the  net  result  is  compound  simple  fragmentation. 
A  computational  approach,  based  on  large-scale  molecular  dynamics,  found  power  laws  with 
n/m  =  — 1 . 15  ±  0.08  ~—l  16  over  a  range  of  one  or  two  orders-of-magnitude  in  fragment  mass. 
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7.3.7  Case  7  {nim  =  5/12) 


Unlike  the  previous  six  cases,  this  final  case  violates  the  law  of  sixths.  Simmons  (1977) 
considered  spray  atomization  in  “a  subset  of  200  tests  selected  almost  at  random  . . .  [from]  over 
2,000  separate  tests  on  about  100  different  nozzle  designs;”  see  also  Spielbauer  et.  al.  (1989).  As 
seen  in  Laney  (2016),  Simmons’  Type  I  root  normal  size  distribution  approximately  corresponds 
to  a  Type  II  Weibull  distribution  with: 


n  5 
m  12 


Ifi 

2l3^2, 


0.41666 


which  is  an  average  of  a  fundamental  and  a  semitone.  Put  another  way,  Simmons  approximately 
obtained  a  quarter  tone.  This  has  long  been  considered  a  universal  size  distribution.  However,  as 
described  below,  the  available  evidence  indicates  that,  by  the  definitions  introduced  here,  it  may 
be  better  described  as  a  compound  size  distribution. 

At  each  fixed  downstream  station,  round  liquid  jets  typically  experience  simple  fragmentation, 
i.e.,  thin  cylindrical  surface  layers  form  ligaments  that  fragment  to  produce  universal  size 
distributions  about  averages  obtained  as  in  Appendix  A.  As  a  result,  spray  atomization  can  and 
often  does  produce  true  universal  size  distributions;  see  also  Laney  (2015a,b,  2016). 

Away  from  the  jet  surface,  liquid  fragments  tend  to  experience  breakup  due  to  aerodynamic 
forces  and  coagulation  due  to  random  collisions;  see,  e.g.,  Ruff  et.  al.  (1992),  Hsiang  &  Faeth 
(1992,  1993),  Chou  et.  al.  (1997).  Such  events  introduce  compounding.  In  addition,  as  described 
by  Laney  (2016),  each  downstream  station  tends  to  fragment  differently.  This  is  because  the  jet 
radius  -  and  thus  the  surface  strain  rate  -  tends  to  increase  linearly  with  distance  from  the  nozzle 
exit.  As  a  result,  averaging  results  across  different  stations  will,  in  many  cases,  yield  compound 
size  distributions.  More  specifically,  when  n  is  different  at  different  stations,  the  average  n  will 
change  continuously  depending  on  the  exact  choice  of  stations  used  to  form  the  average. 

In  Simmons  (1977),  many  if  not  most  of  the  test  results  measured  the  spray  characteristics  away 
from  the  jet  surface,  averaged  across  various  downstream  stations.  This  implies  that 
compounding  occurs  in  many  if  not  most  of  the  individual  tests.  However,  even  if  that  were  not 
the  case,  and  even  if  all  the  individual  tests  involved  simple  fragmentation,  averaging  across 
different  tests  with  different  outcomes  yields  compound  rather  than  universal  size  distributions. 
In  other  words,  when  n  is  different  in  different  tests,  the  average  n  changes  continuously 
depending  on  the  exact  choice  of  tests  used  to  form  the  average. 

While  numerous  papers  have  appeared  in  the  research  literature  providing  additional  evidence 
that  supports  Simmons’  size  distribution,  these  almost  all  involve  averages  across  multiple  tests 
rather  than  specific  tests;  see,  e.g.,  Spielbauer  et.  al.  (1989),  Ruff  et.  al.  (1992),  Hsiang  &  Faeth 
(1992,  1993),  Chou  et.  al.  (1997),  and  Laney  (2015a).  In  other  words,  the  available  evidence 
indicates  that  Simmons  (1977)  obtained  a  compound  size  distribution  rather  than  a  true  universal 
size  distribution.  Like  most  compound  size  distributions,  Simmons’  size  distribution  violates  the 
law  of  sixths. 


32 


7.4  Theory 


Suppose  a  compound  simple  event  can  be  subdivided  into  local  regions  in  time  and  space  where 
the  fragmentation  is  approximately  simple,  i.e.,  where  s  and  m  are  approximately  constant. 
Suppose  the  probability  density  function  for  each  such  local  sub-event  is  as  follows: 

f(X\s\m)  (64) 

Assuming  the  fragments  originating  in  different  regions  in  time  and  space  are  thoroughly  mixed 
together,  and  assuming  m  is  constant,  then  the  probability  density  function  for  the  overall 
fragmentation  event  is  as  follows: 


^w{s)f{X\t,m)ds  (65) 

where  w  is  a  positive  weighting  function  such  that: 

J  w{s)ds  =  1  (66) 

By  the  integral  mean  value  theorem.  Equation  (65)  can  be  written  as  follows: 

^w{s)f{X-,s-,m)ds  =  f{X-,s-,m)  (67) 


where  s  is  an  average  strain  rate.  This  assumes  that  w  >  0  over  the  limits  of  integration  and  that 
/and  w  are  continuous  functions  of  £  ;  see,  e,g..  Theorem  7.2  of  Sahoo  &  Riedel  (1998).  Notice 
that,  because  there  is  no  mean  value  theorem  for  double  integrals,  this  proof  does  not  apply  if  m 
and  £  vary  simultaneously. 

For  example,  if  the  probability  density  function  for  each  local  event  is  a  Type  II  Weibull 
distribution,  then  Equation  (67)  says  that  the  probability  density  function  for  the  overall  mixture 
is  also  a  Type  II  Weibull  distribution  with: 


and: 


(68) 


n  =  n(£,m) 


(69) 


where  n  is  a  either  discrete  or  a  continuous  function  of  £  .  Put  another  way,  compound  simple 
fragmentation  produces  either  universal  or  compound  size  distributions. 
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8.  Piecewise  Simple  and  Piecewise  Compound  Simple  Fragmentation 


8.1  Introduction 


Piecewise  simple  and  piecewise  compound  simple  fragmentation  occurs  when  different  size 
ranges  obtain  different  universal  or  compound  size  distributions.  Section  7.5  found  that 
compound  simple  fragmentation  occurs  when  m  is  constant  and  w  >  0 .  Conversely,  piecewise 
simple  or  piecewise  compound  simple  fragmentation  occurs  when  m  varies  and/or  when  strain 
rate  regions  with  w-0  separate  regions  with  w > 0 . 

As  one  example,  consider  multimode  breakup  of  liquid  droplets;  see,  e.g..  Pilch  &  Erdman 
(1987),  Chou  &  Faeth  (1998),  Dai  &  Faeth  (2001).  In  such  cases,  aerodynamic  forces  cause 
droplets  to  distort  until  they  form  semi-spherical  bags  attached  to  semi-circular  rims;  optionally, 
a  semi-linear  stamen  may  extend  from  the  center  rear  of  the  bag  through  the  center  of  the  rim. 
The  bag,  rim,  and  stamen  probably  each  experiences  simple  fragmentation.  However,  because 
each  feature  experiences  different  strain  rates,  and  the  weighting  function  w  is  zero  except  for 
those  three  discrete  strain  rates,  the  net  result  is  piecewise  simple  fragmentation. 


As  another  example,  consider  explosive  fragmentation  of  thin  cylindrical  shells.  In  many  cases, 
the  largest  fragments  form  in  two-dimensions  while  the  smallest  fragments  form  in  three- 
dimensions.  Then  the  universal  size  distribution  changes  when  the  fragment  dimension  m 
changes.  As  one  specific  case,  suppose  the  universal  changes  from  nlm-113  when  m-3  to 
nlm-112  when  m-2.  These  are  typically  reported  as  a  single  average,  e.g.,  Cohen  (1981) 
found: 


0.433 


1 

2 


(-  -1 

v3^2. 


This  is  about  the  same  as  Simmons  (1977),  as  discussed  earlier  in  Section  7.3.7.  However, 
because  m  is  fixed  at  3  in  Simmons  (1977),  it  occurs  for  different  reason.  As  another  specific 
case,  suppose  the  universal  changes  from  n/ m-2/ 3  when  m-3  to  n/m-l  when  m-2  . 
Again,  these  are  typically  reported  as  a  single  average,  e.g.,  Grady  et.  al.  (2001)  found: 


0.85 


+  1 


This  is  easily  confused  with  the  semitone  n/m-5/6~ 0.83. 

As  a  third  example,  as  discussed  earlier  in  Section  7.6,  Onose  &  Fujiwara  (2004)  studied  porous 
gypsum  spheres  impacted  by  nylon  spheres.  In  each  of  six  reported  tests,  they  found  that  the 
fragment  size  distributions  had  “three  regions  [separated]  by  two  inflections,”  which  they  called 
Region  I,  II,  and  III  for  the  “largest,”  “intermediate”  and  “finest”  sizes,  respectively.  They 
reported  n/ m  for  each  size  region  separately.  In  general,  n/ m  was  substantially  different  in 
different  size  regions.  The  following  subsections  provide  more  examples  of  piecewise  simple 
and  piecewise  compound  simple  fragmentation. 
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8.2  Case  1 


Figure  2  plots  two  Type  II  Weibull  size  distributions  versus  test  data  from  Rosin  &  Rammler 
(1934).  The  smallest  fragments  have  nlm--2l3  while  the  largest  fragments  have  nl  m-213 . 
Brown  &  Woheltz  (1995)  suggest  that  “the  data  consist  of  two  populations:  fines  that 
experienced  a  single  fragmentation  event  and  remained  unaffected  in  spaces  among  larger 
particles  that  were  repeatedly  fragmented  during  milling.” 


Figure  2.  Type  II  Weibull  distributions  vs.  test  data  for  ball-milling  of  iron  from  Rosin  &  Rammler  (1934).  A 
similar  figure  appeared  in  Brown  &  Wohletz  (1995). 


8.3  Case  2 


Wittel  et.al.  (2008)  and  Carmona  et.  al.  (2008)  modeled  collisions  between  three-dimensional 
spheres  using  a  Discrete  Element  Method  (DEM),  where  each  sphere  was  composed  of  a 
hypothetical  material  consisting  of  smaller  “spheres  of  two  different  sizes  . .  .connected  by  beam- 
truss  elements  that  can  elongate,  shear,  bend,  . . .  torque”  and  break.  As  seen  in  Eigure  3,  in  one 
possible  fitting  scheme,  the  largest  fragments  obtained  n/m  «  2  ,  the  middling  fragments 
obtained  n/  mx— 2,  and  the  smallest  fragments  obtained  n!  m^-\. 
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Figure  3.  Type  II  Weibull  size  distributions  vs.  modeling  data  for  colliding  spheres  obtained  by  the  Discrete 
Element  Method  by  Wittel  et.al.  (2008)  and  Carmona  et.  al.  (2008). 

8.4  Case  3 

The  complementary  cumulative  distribution  function  F  is  defined  in  terms  of  the  probability 
density  function /as  follows: 
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F{X)  =  -\f{x)dx-,  f(X)  =  -—  (70) 

This  and  the  following  example  will  use  the  time-proven  technique  of  viewing  F  in  a  plane  in 
which  it  is  linear  or  nearly-linear.  This  view  tends  to  expose  differences  between  size 
distributions  better  than  other  views.  In  Sections  8.2  and  8.3,  the  size  distribution  shifted  from 
positive  to  negative.  Such  dramatic  shifts  are  obvious  in  almost  any  view.  Here  and  in 
Section  8.5,  the  size  distribution  shifts  between  two  positive  size  distributions,  such  that  n/ m 
changes  by  as  little  as  1/12,  which  is  not  apparent  in  many  views. 

o 

To  model  “large-scale  quarry  blastings  of  granitic  gneiss,”  Astrom  et.  al.  (2004)  studied  impact 
fragmentation  of  molded  gypsum  disks  dropped  onto  flat  surfaces  from  heights  of  0.25m  to  10m. 
Each  gypsum  disk  had  a  diameter  of  10.65cm  and  a  thickness  of  2.28cm.  As  shown  in  Figure  4, 
for  a  drop  height  of  0.75m,  the  smallest  fragments  approximately  obey  a  Type  II  Weibull  size 
distribution  (equivalently  a  power  law)  with  nlm--ll3  while  most  of  the  remaining  fragments 
approximately  obey  a  Type  II  Weibull  size  distribution  (equivalently  a  power  law)  with 
nl Similarly,  for  a  drop  height  of  7m,  most  of  the  fragments  approximately  obey  a 
Type  II  Weibull  size  distribution  (equivalently  a  power  law)  with  n/ m--l/3 .  However,  the 
larger  fragments  are  piecewise  compound  simple  or  possibly  completely  non-simple.  Notice  that 
the  fits  suggested  here  differ  from  those  suggested  by  Astrom  et.  al.  (2004). 
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Figure  4.  Type  II  Weibull  size  distributions  vs.  test  data  for  impact  fragmentation  of  molded  gypsum  discs 
dropped  from  heights  of  0.75m  and  7m  taken  from  Astrom  et.  al.  (2004). 


8.5  Case  4 

Mock  &  Holt  (1981,  1983)  gave  highly-detailed  results  for  six  tests  involving  explosively-driven 
cylindrical  metal  shells.  Two  of  these  tests  involved  Armco  iron  shells,  two  of  these  tests 
involved  weak  HF-1  steel  shells,  and  two  of  these  tests  involved  strong  HF-1  steel  shells.  The 
difference  in  the  strengths  of  the  steel  shells  is  due  to  a  difference  in  the  heat  treatment.  The  iron 
shells  had  an  outer  diameter  of  4.5  inches,  an  inner  diameter  of  3  inches,  a  wall  thickness  of  0.75 
inches,  and  a  length  of  8  inches.  The  steel  shells  had  an  outer  diameter  of  4.75  inches,  an  inner 
diameter  of  3  inches,  a  wall  thickness  of  0.875  inches,  and  a  length  of  8  inches.  Notice  that  the 
ratio  of  the  wall  thickness  to  the  outer  radius  is  0.333  for  iron  and  0.368  for  steel,  making  these 
relatively  thick  shells.  The  shells  were  fragmented  by  detonating  2.75kg  cylinder  of  cast-in-place 
Composition  B,  with  57%  inside  the  metal  shells  and  43%  outside  of  the  metal  shells.  The 
charge-to-shell  mass  ratio  was  0.17  for  iron  and  0.15  for  steel. 
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Table  14  summarizes  the  results  of  the  six  tests.  Notice  that  the  curve  fits  given  in  Table  14 
differ  from  those  given  by  Mock  &  Holt  (1981,  1983). 

Table  14.  Summary  of  fragments  size  distributions  vs.  material  strength  for  two  Armco  Iron  and  four  HF-1 


steel  cylindrical  shells  tested  by  Mock  &  Holt  (1981, 1983). 


Exp. 

Casing 

Mat’l 

Casing 

(M 

Strength 

IPa) 

No. 

Size 

Bins 

Type  II  Weibull  Size  Distribution 

Yield 

Tensile 

Smallest  Frags 
{m  =  3) 
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{m  =  2) 
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nim 

n 
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1 
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18 

1/3 

1 

1/2 

1 

1/3 

2 
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It 
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1/3 
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980-1040 
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15 
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3/2 

1/2 

1 
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4 
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It 

65 

1/2 
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5 

Steel 
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67 
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0.58 

6 
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It 

It 

67 

1/2 

3/2 

2/3 

4/3 

0.58 

Figures  5  and  6  show  the  results  of  Experiments  1  and  2  in  a  log-log  plane.  Notice  that  the  two 
experiments  are  identical,  except  for  the  number  of  bins  and  random  shot-to-shot  variation.  In 
Figure  5,  a  Type  II  Weibull  distribution  with  n/m  =  1/3  appears  to  provide  a  good  fit  regardless 
of  fragment  size.  However,  in  Figure  6,  with  eight  times  as  many  bins,  a  Type  II  Weibull 
distribution  n/m  =  1/3  provides  a  good  fit  for  fragments  less  than  600gr,  a  mediocre  fit  for 
fragments  between  600  and  2,000gr,  and  a  poor  fit  for  fragments  greater  than  2,000gr. 


As  an  alternative  to  Figure  6,  Figure  7  shows  the  results  of  Experiment  2  in  a  plane  in  which  F  is 
linear  or  nearly-linear.  In  this  view,  the  size  distribution  clearly  transitions  from  n! m  —  \l 3  for 
the  smallest  fragments  with  m  =  3to  n!  m-H  2  for  the  largest  fragments  with  m-2  .  During 
the  transition,  the  size  distribution  briefly  assumes  the  average: 


n 

m 


1 

2 


{-  -1 


0.41666 


On  the  surface.  Experiments  1  and  2  appear  to  be  textbook  cases  of  simple  fragmentation  with  a 
dimensional  transition.  However,  because  the  iron  shell  is  fairly  thick,  there  is  actually  a 
significant  compounding  effect.  According  to  lower  plot  in  Figure  7,  most  of  the  fragments 
between  about  600  and  2,000gr  are  two-dimensional,  i.e.,  they  retain  two  parallel  smooth 
surfaces  from  the  original  shell.  According  to  upper  plot  in  Figure  7,  these  midrange  two- 
dimensional  fragments  obtain  n/m  =  1/3  just  like  the  smallest  three-dimensional  fragments. 
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Figure  5.  Type  II  Weibull  size  distribution  vs.  test  data  in  a  log-log  plane  for  an  explosively-fragmented  iron 
cylindrical  shell  taken  from  Experiment  I  of  Mock  &  Holt  (1981, 1983).  The  lower  plot  shows  the  fraction  of 
fragments  retaining  portions  of  both  the  inner  and  outer  surface  of  the  original  iron  shell,  which  are  called 
Type  I  by  Mock  &  Holt  (1981, 1983). 


39 


1 


1  10  100  1000  10000 


M{gr) 


Figure  6.  Type  II  Weibull  size  distribution  vs.  test  data  in  a  log-log  plane  for  an  explosively-fragmented  iron 
cylindrical  shell  taken  from  Experiment  2  of  Mock  &  Holt  (1981, 1983).  The  lower  plot  shows  the  fraction  of 
fragments  retaining  portions  of  both  the  inner  and  outer  surface  of  the  original  iron  shell,  which  are  called 
Type  I  by  Mock  &  Holt  (1981, 1983). 
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Figure  7.  Type  II  Weibull  size  distributions  vs.  test  data  in  a  plane  in  wbicb  F  is  linear  or  nearly-linear  for  an 
explosively-fragmented  iron  cylindrical  sbell  taken  from  Experiment  2  of  Mock  &  Holt  (1981, 1983).  The 
lower  plot  shows  the  fraction  of  fragments  retaining  portions  of  both  the  inner  and  outer  surface  of  the 
original  steel  shell,  which  are  called  Type  I  by  Mock  &  Holt  (1981, 1983). 
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Figure  8  shows  the  results  of  Experiment  4  in  a  plane  in  which  F  is  linear  or  nearly-linear.  The 
casing  thickness  is  slightly  larger,  and  the  casing  strength  is  much  larger,  than  in  Experiments  1 
and  2.  As  a  result,  compounding  effects  are  more  pronounced.  In  fact,  the  smallest  fragments 
with  m  =  3  have  n! m  —  \l 2  and  the  largest  fragments  with  m  —  2  have  n/m  =  2/3 ,  reversing 
expectations  based  on  simple  fragmentation.  Notice  that,  during  the  dimensional  transition,  the 
size  distribution  briefly  assumes  the  average: 
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Ifl  2^ 

2  [2^3, 


0.58333 


Eigures  9  shows  the  results  of  Experiment  5  in  a  log-log  plane.  In  this  view,  it  appears  that  a 
Type  II  Weibull  distribution  with  nl  m  —  O.SS  provides  a  good  fit,  except  for  a  few  of  the  largest 
size  bins.  In  other  words.  Experiment  5  appears  to  violate  the  law  of  sixths.  Eor  a  clearer  view, 
Eigure  10  shows  the  results  of  Experiment  5  in  a  plane  in  which  F  is  linear  or  nearly-linear.  In 
this  view,  the  size  distribution  clearly  transitions  from  n/m -1/2  for  the  smallest  fragments  to 
n! m -  2/3  for  the  largest  fragments,  without  assuming  any  intermediate  values  such  as 
n/m  —  0.58 .  In  this  view,  it  is  apparent  that  Experiment  5  satisfies  the  law  of  sixths,  in  a 
piecewise  fashion.  The  material  strength  in  Experiment  5  is  somewhat  less  than  that  in 
Experiment  4.  Comparing  Eigures  8  and  10,  the  main  effect  seems  to  be  eliminating  the 
transition  region  where  n/m  —  0.58 . 

In  general,  a  view  is  a  specific  combination  of  the  following: 

(a.)  a  size  distribution  such  as  f{X),  fyiX),  F{X),or  Fy{X)  for  some  X  and  Y 

(b.)  a  plane  such  as  the  linear-linear  plane,  the  log-log  plane,  or  a  plane  in  which  f(X), 
fy(X),  F(X),or  Fy(X)  is  linear,  piecewise-linear,  or  nearly  so 

(c.)  upper  and  lower  bounds  on  each  axis. 

Comparing  Eigures  6  and  7,  or  Eigures  9  and  10,  shows  how  much  difference  the  view  can  make 
on  the  perceived  goodness  of  fit;  see  also  Eaney  (2015a,  b,  2016).  The  sensitivity  or 
discrimination  of  a  given  view  can  be  determined  by  varying  key  parameters  such  as  n/m.  In 
discriminating  views,  small  changes  in  n/m  lead  to  large  changes  in  the  size  distribution. 
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Figure  8.  Type  II  Weibull  size  distributions  vs.  test  data  in  a  plane  in  which  F  is  Unear  or  nearly-linear  for  an 
explosively-fragmented  steel  cylindrical  shell  taken  from  Experiment  4  of  Mock  &  Holt  (1981, 1983).  The 
lower  plot  shows  the  fraction  of  fragments  retaining  portions  of  both  the  inner  and  outer  surface  of  the 
original  steel  sheU,  which  are  called  Type  I  by  Mock  &  Holt  (1981, 1983). 
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Figure  9.  Type  II  Weibull  size  distribution  vs.  test  data  in  a  log-log  plane  for  an  explosively-fragmented  steel 
cylindrical  shell  taken  from  Experiment  5  of  Mock  &  Holt  (1981, 1983).  The  lower  plot  shows  the  fraction  of 
fragments  retaining  portions  of  both  the  inner  and  outer  surface  of  the  original  steel  shell,  which  are  called 
Type  I  by  Mock  &  Holt  (1981, 1983). 
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Figure  10.  Type  II  WeibuU  size  distributions  vs.  test  data  in  a  plane  in  wbicb  F  is  linear  or  nearly-Unear  for 
an  explosively-fragmented  steel  cylindrical  shell  taken  from  Experiment  5  of  Mock  &  Holt  (1981, 1983).  The 
lower  plot  shows  the  fraction  of  fragments  retaining  portions  of  both  the  inner  and  outer  surface  of  the 
original  steel  shell,  which  are  called  Type  1  by  Mock  &  Holt  (1981, 1983). 


45 


9.  Frequency  Statistics 


To  give  a  sense  for  how  often  various  universal  and  compound  size  distributions  occur,  Table  15 
summarizes  the  results  of  six  fragmentation  studies,  each  of  which  reported  20  or  more  Type  II 
Weibull  size  distributions  (equivalently  power  laws). 


Table  15.  Type  II  Weibull  size  distributions  (equivalently  power  laws)  obtained  in  six  studies  of 
fragmentation  with  »i=3. 


Reference 

Materials 

Type 

■nmin 

■Winax 

No. 

Hartmann  (1969) 

Geological 

Original  Results 

1.80 

3.60 

20 

Turcotte  (1986) 

Geological 

Asteroids 

Meteorites 

Literature  Survey 

1.44 

3.54 

21 

Oddershede  et.  al.  (1993)  (*) 

Gypsum 
Frozen  Soap 
Etc. 

Original  Results 

0.4 

1.9 

20 

Kaminski  &  Jaupart  (1998) 

Magma 

Literature  Survey 

2.9 

3.9 

64 

Onose  &  Fujiwara  (2004) 

Gypsum 

Original  Results 

0.93 

8.46 

28 

Yanovsky,  Tur,  &  Kuklina 
(2010) 

Geological 

Asteroids 

Meteorites 

Glass 

Literature  Survey 

1.44 

4.0 

21 

Total  Number  of  Observed  Fragment  Size  Distril 

buttons 

174 

(*)  Four  size  distributions  were  excluded  because  they  concerned  thin  plates  that  do  not  necessarily  experience 
three-dimensional  fragmentation.  Two  size  distributions  were  excluded  because  they  were  influenced  by  an 
exponential  roll-off  applied  to  the  standard  power  law. 


At  first  glance,  the  exponent  n  in  these  six  studies  seem  to  be  uniformly  random.  However,  upon 
closer  examination,  universals  occur  significantly  more  frequently  than  other  size  distributions; 
see  Figure  11.  More  specifically,  for  \  <—n<  3.5 ,  universals  occur  about  30%  of  the  time. 
Rounding  n  to  the  nearest  0.1,  universals  should  occur  20%  of  the  time  simply  by  random 
chance,  according  to  the  law  of  sixths.  Thus,  for  \  <—n<  3.5  ,  universals  occur  about  50%  more 
often  than  would  be  expected  simply  by  random  chance.  For  -n>  3.5 ,  the  sample  size  is  small, 
and  universals  occur  no  more  often  than  would  be  expected  by  random  chance. 

Due  to  common  curve  fitting  practices,  many  studies  obtain  compound  size  distributions  when 
they  should  have  obtained  universal  size  distributions  and,  vice  versa,  many  studies  obtain 
universal  size  distributions  when  they  should  have  obtained  compound  size  distributions.  The 
most  reliable  studies  are  those  where  conditions  vary,  e.g..  Tables  13  and  14.  If  the  fragment  size 
distributions  change  continuously  with  changing  conditions,  they  are  likely  to  be  compound  size 
distributions.  If  the  fragment  size  distributions  stay  the  same  or  jump  with  changing  conditions, 
they  are  likely  to  be  universals. 


46 


6 


0  0.5  1  1.5  2  2.5  3  3.5  4  4.5  5  5.5  6  6.5  7  7.5  8  8.5  9 


Figure  11.  Exponents  in  Type  11  Weibull  size  distributions  (equivalently  power  laws)  obtained  in  the  six 
studies  listed  in  Table  15.  The  percentages  were  computed  separately  for  each  of  the  six  studies.  The  final 
percentage  was  obtained  by  equally  weighting  the  six  studies  to  obtain  something  akin  to  an  average  PDF. 

Some  questionable  fitting  praetiees  have  already  been  diseussed,  e.g.,  using  a  single  size 
distribution  when  two  or  more  are  aetually  required,  fitting  in  a  single  low-diserimination  view, 
etc.  As  another  example,  researchers  sometimes  assume  the  wrong  kind  of  size  distributions, 
e.g.,  negative  instead  of  positive  size  distributions.  For  example,  Ghosh  et.  al.  (1990)  gave  26 
power  law  indices  for  rock  fragmentation  in  open-pit  copper  mines.  Only  one  is  approximately  a 
universal,  far  less  than  one  would  expect  simply  from  random  chance.  The  likely  source  of  the 
problem  is  seen  in  their  Figure  2,  where  the  chosen  power  law  fit  is  linear  in  the  log-log  plane 
while  the  as-measured  size  distribution  is  highly  non-linear,  i.e.,  the  power  law  is  essentially  a 
tangent  line  to  the  true  size  distribution. 

As  another  example,  closely  related  to  the  previous  one,  the  vast  majority  of  fits  are  done  with 
simple  least  squares  fitting,  rather  than  more  rigorous  statistical  procedures  such  as  Maximum 
Likelihood  Estimation  (MLE),  Pearson’s  chi-square  test,  the  Kolmogorov-Smimov  test,  or  the 
Wald-Wolfowitz  runs  test;  see,  e.g.,  Clauset  et.  al.  (2009),  Gritsevich  et.  al.  (2014). 

As  another  example,  fits  may  be  unduly  influenced  by  measurement  errors  at  the  extremes.  In 
particular,  many  experiments  fail  to  obtain  fully-adequate  statistics  for  the  largest  and  rarest 
fragments,  because  this  may  require  painstaking  and  costly  repetition  of  the  same  test.  Also, 
many  experiments  impose  a  soft  lower  limit  on  fragment  size,  meaning  they  obtain  only  an 
arbitrary  fraction  of  the  fragments  below  a  given  threshold. 

Einally,  fits  may  use  too  many  size  bins,  too  few  size  bins,  fixed  bin  widths,  or  inadequately- 
varying  bin  widths.  Overpopulated  bins  -  which  typically  occur  for  smaller  fragments  -  obscure 
the  characteristics  that  distinguish  one  size  distribution  from  another.  Underpopulated  bins  - 
which  typically  occur  for  larger  fragments  -  tend  to  exaggerate  natural  random  variations.  In  the 
ideal,  bin  widths  should  vary  so  that  all  bins  have  nearly  the  same  number  of  fragments  or,  more 
realistically,  the  number  of  fragments  does  not  change  dramatically  between  adjacent  bins. 
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10.  Conclusions 


This  paper  explains  why  some  studies  have  found  that  different  fragmentation  events  have 
similar  outcomes,  while  other  studies  have  found  that  similar  fragmentation  events  have  different 
outcomes.  More  specifically,  simple  fragmentation  events  -  those  with  simple  geometries  and 
steady,  uniform  strain  rates  -  tend  to  produce  universal  size  distributions.  By  contrast,  complex 
fragmentation  events  -  those  with  complex  geometries  and  unsteady,  non-uniform  strain  rates  - 
tend  not  to  produce  universal  size  distributions.  Having  said  this,  naturally-occurring  mixing  may 
cause  complex  fragmentation  to  resemble  simple  fragmentation,  the  main  difference  being  that 
the  key  size  parameters  may  vary  continuously  rather  than  discretely  with  changing  conditions. 

To  explain  further,  this  paper  finds  that  simple  fragmentation  is  governed  by  Type  II  Weibull 
size  distributions  with  two  free  parameters,  namely, 

n  —  ±1,  ±  2,  ±  3, . . . 
and 

k  =  l,2,3,... 

The  expression  for  n  is  obtained  from  the  three-way  intersection  of  Weibull  size  distributions. 
Gamma  size  distributions,  and  exponential  power  law  size  distributions.  Weibull  and  Gamma 
size  distributions  ensure  geometric  properties  for  fragmentation  and  coagulation,  respectively, 
while  exponential  power  law  size  distributions  ensure  maximum  entropy  properties.  The 
expression  for  k  is  obtained  from  conservation  laws  and  constitutive  relations. 

Simple  fragmentation  is  rarely  seen  in  practice.  Fortunately,  expressions  developed  for  simple 
fragmentation  often  apply  to  non-simple  fragmentation,  after  straightforward  modifications.  The 
following  list  gives  the  five  conditions  required  to  achieve  simple  fragmentation  and,  in  most 
cases,  the  extension  to  non-simple  fragmentation: 

1 .  The  fragment  size  can  be  described  by  a  single  diameter,  such  that  mass  and  surface  area 
are  constant  integer  powers  of  this  diameter.  Counter-examples  include  the  stringy 
highly-distended  fragments  produced  by  non-Newtonian  liquids.  This  assumption  can  be 
relaxed;  see.  e.g.,  the  discussion  of  non-integer  m  in  Section  5  and  non-constant  m  in 
Section  8. 

2.  Material  properties  such  as  density,  surface  tension,  and  fracture  toughness  are 
approximately  constant  throughout  the  fragmentation  event. 

3.  Fragmentation  is  driven  by  a  constant  scalar  tensile  strain  rate.  As  seen  in  Appendix  A.9, 
this  occurs  for  thin  uniformly-radially-expanding  cylindrical  and  spherical  shells, 
uniformly- stretched  bars,  uniformly- stretched  ligaments,  and  so  forth.  This  assumption 
can  be  relaxed.  In  particular.  Sections  7,8,  and  9  show  that  mixtures  resulting  from 
variable  strain  rates,  tensor  strain  rates,  and  non-tensile  strain  rates  often  produce  Type  II 
Weibull  size  distributions  with: 
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n  1112  5, 

m  6  3  2  3  6 


Since  these  are  all  integer  multiples  of  one  sixth,  this  may  be  called  the  law  of  sixths. 

4.  Fragmentation  forces  overwhelm  cohesive  strength,  i.e.,  the  tensile  strain  rate  is 
sufficiently  large.  This  assumption  can  be  relaxed;  see  Section  A.  10. 

5.  Fragmentation  occurs  in  a  single  stage.  In  other  words,  fragmentation  is  not  followed  by 
coagulation  or  additional  fragmentation.  This  assumption  can  be  relaxed;  see,  e.g.,  the 
discussion  of  universal  size  distributions  for  coagulation  in  Section  5  and  the  discussion 
of  multi-stage  fragmentation  in  Section  7.3.5. 

The  non-simple  cases  covered  by  the  various  extensions  listed  above  can  be  categorized  as  either 
compound  simple,  piecewise  simple,  or  piecewise  compound  simple. 

The  last  item  on  the  above  list  merits  further  discussion.  Villermaux  et.  al.  (2004)  argued  that 
liquid  fragmentation  is  often  dominated  by  coagulation.  More  specifically,  they  say  “the 
fragmentation  mechanism  . . .  ,  somewhat  surprisingly,  consists  of  a  coalescence  process  [that]  . . . 
is  in  fact,  generic  of  all  situations  where  drops  come  from  the  capillary  destabilization  of  a 
strongly  corrugated  ligament.”  Their  experimental  results  obtained  a  Type  II  Gamma  size 
distribution  with  b  -3.5  .  As  shown  by  Laney  (2016),  this  approximately  corresponds  to  Type  II 
Weibull  distribution  with  nlm-112  and  m-3,  i.e.,  it  is  universal  just  like  fragmentation  alone 
typically  obtains. 

Similarly,  Michel  et.  al.  (2002)  argue  that  certain  types  of  solid  fragmentation  events  may  also  be 
dominated  by  coagulation.  More  specifically,  they  used  a  Smooth  Particle  Hydrodynamics  (SPH) 
code  to  simulate  fragmentation  of  monolithic  basalt  bodies,  meant  to  represent  large  asteroids, 
and  an  N-body  code  to  simulate  coagulation.  They  conclude  “that  the  parent  body  is  first 
completely  shattered  at  the  end  of  the  fragmentation  phase,  and  then  subsequent  gravitational 
reaccumulations  lead  to  the  formation  of  an  entire  family  of  large  and  small  objects ....”  They 
obtained  ten  power  laws  (equivalently  Type  II  Weibull  size  distributions  with  n  <  0),  three  of 
which  are  approximately  fundamentals  and  four  of  which  are  approximately  semitones,  much 
greater  than  would  be  expected  by  random  chance.  Again,  the  combination  of  fragmentation  and 
coagulation  obtains  the  same  universal  size  distributions  as  fragmentation  alone. 

Many  researchers  have  stressed  the  role  of  secondary,  tertiary,  etc.  fragmentation,  in  which  the 
original  ‘parent’  fragments  break  up  into  smaller  ‘child’  fragments  due  to  collisions, 
aerodynamic  forces,  and  so  forth.  Similarly,  there  may  be  multiple  stages  of  coagulation.  This 
treatment  helps  to  explain  why  multiple  stages  of  fragmentation  and  coagulation  often  produce 
universal  size  distributions  just  like  simple  one-stage  fragmentation. 

In  short,  this  paper  finds  that  complex  multistage  events  -  which  produce  a  diverse  mixture  of 
fragments  formed  at  different  times,  in  different  places,  in  different  dimensions,  and  under 
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different  conditions  -  may  often  be  treated,  in  aggregate,  as  a  single  virtual  event  governed  by  a 
single  universal  or  compound  size  distribution. 

Where  universal  size  distributions  apply,  a  key  unsolved  problem  is  choosing  n  and  k.  One 
common  approach  is  to  assume  frequently-observed  values  such  as  n  =  1  or  2  and  k  -  2  or  3.  A 
better  but  less  common  approach  is  to  estimate  n  and  k  from  experimental  databases  or  curve  fits 
to  experimental  results,  e.g.  Cunningham  (2005).  A  third  possible  approach  is  stability  analysis, 
as  discussed  briefly  in  Section  7.1.  This  is  common  when  easily-visible  oscillations  occur  prior 
to  fragmentation  but  may  apply  more  generally. 
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Appendix  A.  Algebraic  Conservation  Laws  for  Simple  Fragmentation  Events 
A.l  Introduction 

This  appendix  derives  simple  algebraic  expressions  for  average  fragment  sizes  from  conservation 
of  mass,  momentum,  and  energy,  plus  constitutive  relations.  This  treatment  assumes  that 
fragmentation  is  driven  by  a  scalar  tensile  component  of  the  strain  rate  or,  alternatively,  the 
tensile  component  of  the  kinetic  energy.  This  treatment  builds  on  earlier  proofs  including  the 
following: 

Spherical  Geometries  Grady  &  Kipp  (1993) 

Cylindrical  Geometries  Mott  &  Linfoot  (1943);  Grady  (1982);  Wu  et.  al.  (1992) 
Cartesian  Geometries  Lefebvre  (1992a,b);  Kirane  et.  al.  (2015) 

The  current  approach  has  several  advantages  over  earlier  approaches.  First,  earlier  approaches 
usually  did  not  commit  to  a  particular  average  fragment  size,  such  as  the  mass  mean  diameter; 
the  current  approach  is  specific.  Second,  earlier  approaches  usually  assumed  that  the  original 
fragmenting  body  was  a  solid  or  hollow  cylinder,  sphere,  cube,  and  so  forth;  this  approach  makes 
no  such  assumption.  Third,  earlier  approaches  usually  assumed  that  the  fragments  were  spheres, 
cubes,  and  so  forth;  this  approach  makes  no  such  assumptions.  Finally,  earlier  approaches 
assumed  a  constant  multiplicative  factor;  this  approach  shows  that  the  multiplicative  factor  is  not 
constant  but  rather  depends  on  the  fragment  size  distribution,  fragment  shape,  the  choice  of 
average  fragment  diameter,  and  the  minimum  fragment  size  inherent  in  most  experimental 
measurements. 

A.2  Algebraic  Forms  of  Conservation  of  Energy 

If  the  two  primary  forms  of  energy  are  kinetic  energy  and  surface  energy,  conservation  of  energy 
requires: 


(A.l) 


where  V  is  the  tensile  component  of  velocity  before  fragmentation  measured  in  units  such  as 
cm/s,  f/,  is  the  tensile  component  of  velocity  of  fragment  i  measured  in  units  such  as  cm/s, 
is  the  surface  area  before  fragmentation  measured  in  units  such  as  cm^,  is  the  surface  area  of 
fragment  i  measured  in  units  such  as  cm  ,  M  ■  is  the  mass  of  fragment  i  measured  in  units  such 
as  g,  F;  is  the  surface  energy  per  unit  area  of  fragment  i  measured  in  units  such  as  dyne/cm,  and 
N  is  the  total  number  of  fragments.  Notice  that: 
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(A.2) 


III  '-V^-dM^±^VfM, 

body  j—^ 

In  other  words,  the  tensile  kinetic  energy  of  the  whole  is  equal  to  the  tensile  kinetic  energy  of  its 
parts.  Also  assume  that: 


(A3) 

1=1 

In  other  words,  assume  the  surface  energy  of  the  whole  is  much  less  than  the  surface  energy  of 
its  parts.  Finally  assume  that: 


Uj  «V.  (A.4) 

In  other  words,  assume  that  the  vast  majority  of  tensile  kinetic  energy  is  converted  to  surface 
energy.  Notice  that  Equation  (A.4)  concerns  only  the  tensile  component  of  velocity.  Fragments 
generally  retain  a  substantial  non-tensile  component  of  velocity;  see,  e.g.,  Kennedy  (1998). 

With  the  above  assumptions.  Equation  (A.2)  can  be  written  as  follows 

(A.5) 

^  <=1 

Equation  (A. 5)  allows  simple  analytical  solutions  if  V ^  is  a  power  of  Dj  and  E,.  is  a  power  of 
D. .  More  specifically,  suppose: 


=  Vo 

(A.6) 

r,  -  r„  D! 

(A.7) 

where  Vq,  v  ,  Eq  and  y  are  constants.  In  addition.  Equation  (1)  says: 

^,=c.-,or'  (A.8) 

(A.9) 

Notice  that  if  m  =  1 ,  then  c„_i  =  Cq  =  2<5^ .  In  other  words,  <5^  =  2^^  where  <%o  is  a  fixed  cross- 
sectional  area.  Equation  (A. 5)  becomes: 
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(A.  10) 


N  *  9^  F 

Zr^lv-y+l  j _ _  0 

^  ~  r  nV^ 

i=l  J^r+m-l  ^mP^Q 


By  Equation  (5),  this  becomes: 


2c^-iro 

Cmpyo 


2f  —  7  +  l>0,  7  +  m  — 1>0 


(A.  11  a) 


y'  _  ^mP^O 

9^  F 

0 


2f  —  7  +  l<0,  7  +  m  — 1>0 


(A.  11b) 


where: 


A  =  D 


(A12) 


Y  =  D 


(A.13) 


If  7  +  m-l  =  0  then  is  replaced  by  and  Xy^^^  is  replaced  by  A'^^. 

This  is  the  final  result.  Notice  that  Equation  (A.  11a)  applies  primarily  to  positive  size 
distributions,  e.g.,  Type  II  Weibull  distributions  with  n>0.  Similarly,  Equation  (A.  1  lb)  applies 
primarily  to  negative  size  distributions,  e.g..  Type  II  Weibull  distributions  with  n  <0 . 

A.3  Case  1 

Starting  with  the  easiest  case,  suppose  that: 


Vj  =  Vq  =  const. 


(A.  14) 


r,  =  Eq  =  const. 


(A.15) 


Using  conservation  of  momentum.  Equation  (A.  14)  can  be  written  as  follows: 


Vj=V,  =  -eD, 


(A.16) 


where  is  a  fixed  geometric  length  scale  measured  in  units  such  as  cm  and  £  is  a  tensile 
strain  rate  measured  in  units  such  as  1/s.  Eor  example,  Grady  (1981)  suggests  0^=  L  where  L  is 
the  initial  length  of  a  linear  body.  Other  possibilities  include:  the  thickness  of  a  cylindrical  shell, 
e.g.,  Mott  et.  al.  (1944);  the  thickness  of  a  flat  sheet,  e.g.,  Eefebvre  (1992a,  b);  and  the 
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characteristic  size  of  or  spacing  between  pre-existing  or  strain-induced  features  such  as  flaws, 
cracks,  voids,  bubbles,  or  instability-driven  oscillations,  e.g.,  Curran  &  Seaman  (1996),  Zhou  et. 
al.  (2005). 

Using  constitutive  equations.  Equation  (A.  15)  can  be  written  as  follows: 

r,=a  (A.  17) 


for  in  viscid  liquids  where  cr  is  the  surface  tension  measured  in  units  such  as  dyne/cm;  see,  e.g., 
Grady  (1982)  and  Grady  &  Kipp  (1982).  Similarly: 


(A.18) 


for  solids  where  is  the  work  of  fracture  measured  in  units  such  as  dyne/cm.  Zhou  et.  al. 

(2005)  describe  as  the  “fracture  energy  dissipated  by  the  crack.”  Notice  that: 

g,=kIie 


E-pc] 


1/2  2 

where  is  the  fracture  toughness  measured  in  units  such  as  g/(cm  s  ),  E  is  the  Young’s 
modulus  measured  in  units  such  as  dyne/cm  ,  and  c,  is  the  elastic  wave  velocity  measured  in 
units  such  as  cm/s.  Then  Equation  (A.18)  can  be  rewritten  as  follows: 


(A.  19) 


Substituting  V^-  s 0^12,  l>  =  0  and  y  =  0  into  Equations  (A.  11a),  (A.  12),  and  (A.  13)  yields: 


D, 


c^pDls^ 


(A.20) 


Suppose  that  Equation  (A.  14)  is  replaced  by  the  following: 

r 

V,  =  const. - - —  (A.21) 

^  psD,D^ 

while  Equation  (A.  15)  is  the  same  as  before.  Substituting  Vq  =  const.  Eq  /{peD^) ,  v  =  -\  and 
y  =  0  into  Equations  (A.  11b),  (A.  12),  and  (A.  13)  yields: 
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(A.22) 


D' 


^  avg 


const. 


Notice  that  Equation  (A. 20)  applies  to  positive  size  distributions  while  Equation  (A.22)  applies 
to  negative  size  distributions. 

Equations  (A. 20)  and  (A.22)  can  be  converted  to  other  average  sizes  using  Equations  (1 1)  to 
(14).  In  general,  this  produces  expressions  for  average  fragment  diameters  in  the  following  form: 

const. — (A.23) 

pDle^ 


Eor  example,  for  positive  size  distributions  for  liquids,  substituting  Equation  (A.  16)  and  (A.  17) 
into  Equation  (A.23)  gives: 


const.  (A. 24) 

which  approximately  agrees  with  well-known  results  for  atomization  of  flat  and  conical  liquid 
sheets;  see,  e.g.,  Eefebvre  (1992a,  b)  and  Barreras  et.  al.  (2006). 

A.4  Case  2 

Suppose  that: 


V^P-sD,  (A,25) 

Notice  that  Equation  (A. 25)  is  the  same  as  Equation  (A.  16),  after  replacing  by  D . ,  In 

addition,  based  on  constitutive  equations  for  solids,  Grady  (1988)  and  Grady  &  Kipp  (1993) 
suggest: 


(A.26) 

D 

where  7  is  a  constant  critical  failure  strength,  such  as  the  tensile  yield  strength,  measured  in 
units  such  as  dynes/cm  and  is  a  unitless  “critical  volumetric  strain  at  void  coalescence” 
assumed  to  be  0.15.  Substituting  Vg  =  f/2,  =  1,  Eg  =  s^Y  /  6  ,  andy  =  1  into  Equations 

(A.  11a),  (A.  12),  and  (A.  13)  yields: 


(£>") 


D" 


avg 


(A.27) 
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Suppose  Equation  (A. 25)  is  replaced  by  the  following: 


r,. 

V.  =  const. — ^  (A.28) 

while  Equation  (A. 26)  is  the  same  as  before.  Substituting  Equation  (A. 26)  into  Equation  (A.28) 
yields: 


Vj  =  const. 


(A.29) 


Substituting  Vq  -  const.  s^Y l{ps) ,  u  =  -\,  Eq  =  / 6  ,  and y  =  \  into  Equations  (A.  1  lb), 

(A.  12),  and  (A.  13)  yields: 


(£>')'  =  const. 


(A.30) 


Notice  that  Equation  (A. 27)  applies  to  positive  size  distributions  while  Equation  (A.30)  applies 
to  negative  size  distributions. 

Equations  (A. 27)  and  (A.30)  can  be  converted  to  other  average  sizes  using  Equations  (1 1)  to 
(14).  In  general,  this  produces  expressions  for  average  fragment  diameters  in  the  following  form: 


Z  \  1/2 

const.  sY  ' 


ps^ 


(A.31) 


Eor  example,  for  positive  size  distributions  with  m  =  3  ,  Equation  can  be  written  as 

follows: 


‘'M  avg 


or: 


where: 


^avg  - 


•2 


pa  ) 


C 


4^2 


(A.32) 


(A.33) 


(A.34) 
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For  spherical  fragments: 


^  ^  =  6  (A.35) 

C3  n  16 

Then: 

C  =  — ^  (A.36) 

If  «  1 ,  this  agrees  with  expressions  given  by  Mott  (1947),  Grady  (1982),  Grady  (1988), 

and  Grady  &  Kipp  (1993)  for  explosively-driven  metal  shells. 

A.5  Case  3 


Suppose  that: 


V  =-iD. 

j  2  ■' 

(A.37) 

II 

II 

8 

s 

(A.38) 

Notice  that  Equation  (A. 37)  is  the  same  as  Equation  (A. 25)  and  Equation  (A. 38)  is  the  same  as 
Equation  (A.  15).  Substituting  Vq  -  sI2,  u  =  1,  and  /  =  0  into  Equations  (A.  1  la),  (A.  12),  and 

(A.  13)  yields: 

/r)3\  _  ^^m-1^0 

(A.39) 

Suppose  Equation  (A. 37)  is  replaced  by  the  following: 

r,. 

V,  =  const. — — 
peD] 

(A.40) 

while  Equation  (A. 38)  is  the  same  as  before.  Notice  that  Equation  (A. 40)  is  the  same  as  Equation 
(A. 28).  Substituting  Vg  -  const.  Eg  l{sp) ,  v--2,  and  y  =  0  into  Equations  (A.  11b),  (A.  12), 

and  (A.  13)  yields: 

(D^)'  =  const. 

(A.41) 
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Notice  that  Equation  (A. 39)  applies  to  positive  size  distributions  while  Equation  (A.41)  applies 
to  negative  size  distributions. 

Equations  (A. 39)  and  (A.41)  can  be  converted  to  other  average  sizes  using  Equations  (1 1)  to 
(14).  In  general,  this  produces  expressions  for  average  fragment  diameters  in  the  following  form: 


/  ^  \l/3 

constr„  ' 


pe 


(A.42) 


which  agrees  with  the  well-known  expressions  given  by  Mott  &  Einfoot  (1943),  Grady  (1982), 
and  Grady  &  Kipp  (1993)  for  explosively-driven  cylindrical  metal  shells. 

Eor  example,  for  positive  size  distributions  with  m  =  3  ,  Equation  (A.  39)  can  be  written  as 
follows: 


M 


SCjEo 


^  avg 


(A.43) 


Because  this  example  is  so  well-known,  more  details  will  be  given  for  this  example  than  for  the 
others.  By  Equations  (11)  to  (14),  Equation  (A.43)  can  be  written  as  follows: 


/  \  1/3 

^  CT  ^ 


ps^  j 


(A.44) 


D'  = 

avg 


ps  J 


where: 


(A.45) 


C  = 


8c^ 


C'  = 


8^2 


(A.46) 


(A.47) 


Eor  Type  II  Weibull  size  distributions.  Tables  3  and  4  give: 


A 


Qpi  ~ 


i  +  - 
V  nj 


(.  A 

1  +  - 

E 

1+^ 

V  nj 

^  nJ 
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f 


R  =  r 


i  +  - 
V  nj 


f 


1-- 
V  n) 


Then: 


f 


- 


1+- 
V  n) 


f 


1 


1  + 

V  n) 


8c2  r| 


c  = 


1+ ^ 

V  n) 


r  2^' 


1  + 

V  n) 


C3r| 


i  +  - 
V  n) 


c  = 


3 

f  23 

3 

r 

r 

(  nj 

^  nj 

CjTl 


1- ^ 

V  nj 


1  +  - 

V  nJ 


Similarly,  by  Equations  (11)  to  (14),  Equation  (A. 43)  can  be  written  as  follows: 


D 


X  avg 


) 


^Xavg  ~ 


C'T 

) 


Then: 


^  %c,Q\ 

^  c,Q%S^ 


C  = 


Eor  Type  II  Weibull  size  distributions.  Tables  3  and  4  give: 


(A.48) 

(A.49) 

(A.50) 

(A.51) 

(A.52) 

(A.53) 
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Then: 


Qpi  ~ 


r  1  + 


\\  2 

r  i+-  r  i  +  - 

\  n  \  n 


Qm  = 


r  1+ 


f  1 W  3 

r  i  +  -  r  i+- 


f  3W  1 

r  i  +  -  r  i+- 

\  nj  {  n 


r  1+ 


A\  2 

r  i+-  r  i  +  - 


r  1+ 


(  i+3l  f  i 

r  i+—  r  i+- 


- 


r  1+ 


8c2r  1 + 


= 


c,r  i  +  -  r  1+ 


8c2r  1  + 


c  - 

'-'5?  “ 


c,r  i  +  -  r  1+ 


(A.54) 


(A.55) 
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Sc^r 


c  = 


(  2 
1  +  - 

V  nj 


1  +  - 
V  nj 


f 


CjF 


lA 


1  +  - 
V  n) 


f 


1  +  - 
V  n) 


8c2r  r| 


c  = 


1+- 
V  nj 


C3r| 


2 

1  +  - 

r 

1  +  - 

V  nj 

^  n) 

(A.56) 


(A.57) 


Notice  that: 


8co 


=  8 


(A.58) 


for  rectangular  fragments  and: 


8C2  _  8;t 
C3  n  !  6 


(A.59) 


for  spherical  fragments.  For  spherical  fragments,  Figure  A.l  shows  how  C,  ,  and  vary 

with  the  Weibull  exponent  n,  where  small  n  corresponds  to  large  fragment  size  spreads  and  large 
n  corresponds  to  small  fragment  size  spreads. 
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n 


Figure  A.l.  Multiplicative  constant  vs.  exponent  n  in  a  Type  II  Weibull  size  distribution  for  spherical 
fragments  where  “number  weighted”  refers  to  Equation  (A.48),  “surface  area  weighted”  refers  to  Equation 
(A.54),  and  “mass  weighted”  refers  to  Equation  (A.56). 

As  seen  in  Figure  A.l,  the  number- weighted  average  fragment  diameter  changes  dramatically 
with  n.  More  specifically,  C  increases  by  a  factor  of  over  150  when  n  increases  from  1  to  oo. 
This  is  mainly  because,  for  large  fragment  size  spreads,  number- weighted  averages  are  heavily 
skewed  by  a  very  large  number  of  very  small  fragments.  In  fact,  the  very  small  fragments  may 
contain  a  significant  fraction  of  the  overall  surface  energy,  even  though  they  contain  only  a 
trivial  fraction  of  the  overall  mass,  momentum,  and  kinetic  energy. 

As  seen  in  Figure  A.l,  surface-area-  and  mass-weighted  averages  do  not  change  dramatically 
with  n.  More  specifically,  increases  by  a  factor  of  about  2.2  and  decreases  by  a  factor  of 
about  0.8  when  n  increases  from  1  to  oo.  This  is  mainly  because,  for  large  fragment  size 
spreads,  surface-area  and  mass-weighted  averages  are  skewed  -  albeit  not  dramatically  so  -  by  a 
very  small  number  of  very  large  fragments. 

Typical  experiments  obtain  and  C  where  the  overbar  refers  to  the  effects  of  a  minimum 
fragment  size  .  In  cases  where  is  unusually  small,  number- weighting  with  a  minimum 
fragment  size  may  roughly  correspond  to  surface-area  weighting.  In  other  words: 
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For  example,  for  spherical  fragments  with  Type  II  Weibull  size  distributions: 


C  =21.6  for  n  =  I 


(A.60) 


If  is  not  unusually  small,  C  may  roughly  equal  or  something  larger. 


Previous  treatments  typically  assumed  C  is  constant.  The  fact  that  C  is  not  constant,  but  rather 
depends  on  the  fragment  size  distribution  through  n,  on  the  fragment  shape  through  m  and 
Cm^i! Cm  ,  and  on  the  minimum  fragment  size,  helps  to  explain  the  wide  range  of  values  for  C 

found  in  the  research  literature.  For  example,  Grady  (1982)  suggests  C  =  20  while  Grady  & 
Kipp  (1993)  suggest  C  =  24,  48,  or  720. 

A.6  Case  4 


Suppose  that: 


8D] 

V.  =  ' 

^  2Do 

(A.61) 

r  _ 

(A.62) 

Notice  that  Equation  (A. 61)  is  the  same  as  Equations  (A. 25)  times  Dj  / Dq.  Similarly,  Equation 
(A. 62)  is  the  same  as  Equation  (A.  19)  times  D.  / D^.  Substituting  =  s u  =  2, 

Fq  =  I{2pc]Dq)  ,  and  y  =  \  into  Equations  (A.  11a),  (A.  12),  and  (A.  13)  yields: 


^  avg 


CmP 


Suppose  Equation  (A. 61)  is  replaced  by  the  following: 


V  ■  =  const. 


peD] 


(A.63) 


(A.64) 


while  Equation  (A. 62)  is  the  same  as  before.  Notice  that  Equation  (A.64)  is  the  same  as  Equation 
(A. 28)  times  .  Substituting  Equation  (A. 62)  into  (A.64)  yields: 


V  ■  =  const. 


^2  2 

P  c 


sD] 


(A.65) 
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Substituting  Vq  -  const.  l{p^c]s) ,  v  =  —2,  Tq=  KlpclDo) ,  and  y  =  \  into  Equations 
(A.  11b),  (A.  12),  and  (A.  13)  yields: 


(D'*)'  =  const.- 

^  ^  D  avg 


CmKtD, 


(A.66) 


Notice  that  Equation  (A. 63)  applies  to  positive  size  distributions  while  Equation  (A.66)  applies 
to  negative  size  distributions. 

Equations  (A. 63)  and  (A.66)  can  be  converted  to  other  average  sizes  using  Equations  (11)  to 
(14).  In  general,  this  produces  expressions  for  average  fragment  diameters  in  the  following  form: 


fconst..S:^7^^  ' 


PCs^ 


J 


(A.67) 


which  approximately  agrees  with  Equation  (13.25)  in  Curran  &  Seaman  (1996)  for  explosively- 
driven  expanding  cylindrical  metal  shells,  where  is  the  “characteristic  size  of  the  microcrack 

distribution  . . .  [for]  flaws  (cracks)  that  have  the  potential  of  growing  to  coalescence  to  form 
fragments.” 

A.7  Case  5 


Suppose  that: 


= 

(A.68) 

2D, 

r,  =  Eg  =  const. 

(A.69) 

Notice  that  Equation  (A. 68)  is  the  same  as  Equation  (A. 61)  while  Equation  (A.69)  is  the  same  as 
Equation  (A.  15).  Substituting  -  s I{2Dq)  ,  v  =  2,  and  y  =  0  into  Equations  (A.  11),  (A.  12), 
and  (A.  13)  yields: 


(D^) 


Suppose  Equation  (A.69)  is  replaced  by  the  following: 


EnOn 

y,  =  const.^^ 


psD. 


(A.70) 


(A.71) 
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while  Equation  (A. 69)  is  the  same  as  before.  Notice  that  Equation  (A.71)  is  the  same  as  Equation 
(A. 64).  Substituting  Eg  =  const. E^Dq /(/:»£) ,  u  =  -3,  and  y  =  0  into  Equations  (A. 11),  (A. 12), 
and  (A.  13)  yields: 


(D^)'  =  const. 

D  avg 


(A.72) 


Notice  that  Equation  (A. 70)  applies  to  positive  size  distributions  while  Equation  (A.72)  applies 
to  negative  size  distributions. 

Equations  (A.70)  and  (A.72)  can  be  converted  to  other  average  sizes  using  Equations  (11)  to 
(14).  In  general,  this  produces  expressions  for  average  fragment  diameters  in  the  following  form: 


const  EqDq  ^ 

V  . 


1/5 


(A.73) 


Eor  example,  by  Equation  (A.  18),  suppose  that: 


In  addition,  suppose  that: 


D 


0 


(A.74) 


(A.75) 


Zhou  et.  al.  (2005)  call  this  the  “critical  opening  distance  of  the  nucleated  crack.”  Substituting 
Equations  (A.74),  (A.75),  and  p  =  E/c^into  Equation  (A.73)  yields: 


^  const 
E^'^Ys 


(A.76) 


which  agrees  to  some  extent  with  Equation  (12)  in  Zhou  et.  al.  (2005)  for  dynamic  fragmentation 
of  a  one-dimensional  brittle  bar;  see  also  Drugan  (2001).  However,  instead  of  a  strain  rate 
exponent  of  —  2  /  5  =  —0.4 ,  Zhou  et.  al.  (2005)  obtained  —  0.4264. 

A.8  Case  6 


Suppose  that: 


sD] 

V:  =  '' 


2Dl 


(A.77) 
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r,  = 


(A.78) 


G.D, 

2D, 

Substituting  =  ^/(2Dq  ),  u  =  5,  T,  -  K'^D,)  and  /  =  l  into  Equations  (A.  11),  (A.  12),  and 

(A.  13)  yields: 


D"‘  avg 


(A.79) 


In  general,  this  produces  expressions  for  average  fragment  diameters  in  the  following  form: 


^  const  .G^Dq  ^ 

ps^  j 


(A.80) 


and  expressions  for  average  fragment  masses  in  the  following  form: 


(A.81) 


which  agrees  to  some  extent  with  Equation  (6.7)  in  Mott  et.  al.  (1944)  for  explosively-driven 
expanding  cylindrical  metal  shells.  However,  instead  of  a  strain  rate  exponent  of 
—  2x3/10  = —0.6,  Mottet.  al.  (1944)  obtained  —0.623. 

A.9  Tensile  Strain  Rates 


The  examples  given  above  require  knowing  a  scalar  tensile  strain  rate,  which  is  assumed  to  be 
large,  steady,  and  uniform.  The  best  known  cases  involve  thin  rapidly-expanding  rings  and 
cylindrical  shells  of  radius  ^  where: 


8 


(A.82) 


according  to  Hoggatt  &  Recht  (1968);  see  also  Warnes  et.  al.  (1985)  and  Jones  et.  al.  (2013).  Eor 
thin-rapidly  expanding  spheres  of  radius  : 


(A.83) 


according  to  Al-Hassani  &  Johnson  (1969);  see  also  Elorca  &  Juanicotena  (1997)  and  Buy  & 
Elorca  (2002). 
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Notice  that  the  previous  two  expressions  apply  only  to  thin  fragmenting  objects.  At  a  radius  r  in 
a  thick  rapidly-expanding  ring  or  cylindrical  shell  with  an  outer  radius  ^  : 


s  = 


where  =  rr  so  that: 

r 

s  =  - 
r 


(A.84) 


(A.85) 


Similarly,  at  a  radius  r  in  a  thick  rapidly-expanding  spherical  shell  with  an  outer  radius  ^  : 


£  =  2 


where  =  rr^  so  that: 


£ 


2r 

r 


(A.86) 


(A.87) 


Finally,  for  uniformly-stretched  one-dimensional  objects  such  as  rods,  bars,  strings,  filaments, 
threads,  ligaments,  etc.  of  length  L: 


.  L 

£-  — 
L 


A.IO  Extensions  to  Low  Strain  Rates 


(A.88) 


Many  fragmentation  events  violate  Equation  (A. 3).  For  example,  Edwards  &  Deal  (2011)  found 
explosively-driven  aluminum  alloy  cylindrical  shells  that  produced  as  a  few  three  fragments.  For 
another  example,  in  a  literature  survey,  Zhou  et.  al.  (2006c)  found  expanding  ring  experiments 
that  produced  as  few  as  18  fragments.  The  expressions  given  in  Section  A. 2  can  be  extended  to 
lower  strain  rates  by  eliminating  Equation  (A. 3),  i.e.,  by  accounting  for  the  surface  area  of  the 
original  fragmenting  body. 

First,  consider  thin  fragmenting  objects  such  as  shells  and  sheets.  Using  conservation  of  mass, 
Eefebvre  (1992a)  argues  that: 


(A.89) 
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where  the  shell  or  sheet  thickness  t  is  assumed  to  be  small.  The  factor  of  two  accounts  for  the 
two  sides  of  the  shell  or  sheet.  In  this  case,  Equation  (A.  10)  becomes: 


N 

Z-l  i  N  '  .^t/2 

j=l  j=l 


AT  ^  j^r+nt  1  Y 

N  ~  ^  Tr2 


1°.  7+m— 1 


i=\ 


/=! 


(A.90) 


Applying  Equation  (5),  Equation  (A.  11)  becomes: 


4r  2c  Y 
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where: 
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Equation  (A.91)  has  a  simple  solution  if: 
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In  this  case: 
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If  Vq  =  sDq  /  2 ,  then  Equation  (A.93a)  is  the  same  as  Equation  (A. 20)  for  large  enough  strain 
rates.  If  /  =  0,  r-EQ-cr,  m-3,c^=7r/6,  and  =  tt  .  Then  Equation  (A.93a)  becomes: 


D  avg 

t 


4cr 


(A.94) 


which  is  similar  to  Eefebvre  ( 1992a, b). 

Alternatively,  for  uniformly- stretched  one-dimensional  objects  such  as  expanding  rods,  bars, 
strings,  filaments,  threads,  ligaments,  etc.  of  initial  length  and  fixed  cross-sectional  area 
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assuming  /  =  0  and  m  =  1 .  The  same  expression  applies  to  expanding  rings  if  =  2?!^^ . 
Equation  (A.95)  assumes  that: 
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(A.96) 


Eor  example,  as  in  Section  A. 3,  suppose  that  =  0 .  Then  Equation  (A.95)  becomes: 


1  c, 
— -r- 


.pv^ 


4r„, 


(A.97) 


Note  that  Equation  (A. 99)  is  formally  similar  to  Equations  (A.93)  and  (A.94).  As  another 
example,  as  in  Section  A. 5,  suppose  that  Vq  -  sI2  and  v  =  \.  Then  Equation  (A.95)  becomes: 
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By  the  cubic  equation: 
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Eor  large  strain  rates,  this  is  the  same  as  Equation  (A. 39)  with  m  =  1  and  c^_i  =Cq=  2^  . 
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In  addition  to  the  approach  described  here,  Drugan  (2001),  Zhou  et.  al.  (2004,  2005,  2006a,b,c), 
Levy  &  Molinari  (2010)  and  Levy  et.  al.  (2012)  compare  several  other  approaches  for  modeling 
low-strain-rate  fragmentation  with  m  =  1 .  For  example.  Levy  &  Molinari  (2010)  suggest  an 
approximate  solution,  formally  similar  to  Equations  (A. 95)  and  (A.97),  found  using  empirical 
curve  fits  to  models  and  test  data.  The  comparisons  given  in  these  references  help  show  how 
large  the  strain  rate  needs  to  be  before  Equation  (A.  1 1)  applies. 

A.ll  Conclusions 

This  treatment  assumes  that  a  dominant  scalar  tensile  strain  rate  s  is  steady  and  uniform.  While 
this  is  a  severe  restriction,  many  problems  can  be  subdivided  into  local  regions  where  the  tensile 
strain  rate  is  approximately  steady  and  uniform.  Eor  example,  as  seen  in  Section  A. 9,  the  strain 
rate  is  steady  and  uniform  in  each  radial  layer  in  a  thick  rapidly-expanding  solid  ring,  cylinder,  or 
sphere.  More  complex  cases  require  applying  expression  such  as  Equation  (A.ll),  (A.91),  or 
(A. 95)  locally  in  first-principles  structural,  fluid  dynamics,  or  hydrodynamics  codes;  see,  e.g., 
Wilson  et.  al.  (2001). 

If  the  tensile  strain  rate  £  is  large  -  in  addition  to  being  steady  and  uniform  -  then  a  series  of 
examples  shows  that: 


(A.lOl) 

where  k  is  any  positive  integer  and  the  subscript  ref  refers  to  the  averages  described  in  Section  3. 
In  addition,  the  average  fragment  diameter  depends  on  the  material  type  (e.g.,  solid  vs.  liquid), 
the  material  properties  (e.g.  density,  strength,  elasticity),  the  fragment  size  distribution,  the 
fragment  shape,  and  the  minimum  fragment  size,  if  any.  The  resulting  expressions  are  well- 
validated  for  k  equal  to  1,2,  and  3  but  are  less  well-validated  for  larger  k.  Despite  the  fact  that 
the  earliest  work  by  Sir  Nevill  Mott  and  his  co-workers  found  a  variety  of  different  values  for  k, 
the  modern  research  literature  has  focused  predominantly  on  k  -  3. 

The  approaches  described  here  are  perhaps  more  popular  in  theory  than  in  practice.  This  is  partly 
because  they  may  predict  average  fragments  sizes  near  or  even  below  the  minimum  fragment 
sizes  imposed  in  typical  experiments;  this  is  especially  true  for  the  count  mean  diameter. 
Common  semi-empirical  alternatives  include:  the  Gumey-Sarmousakis  equation  for  cased 
munitions,  e.g.,  Victor  (1996)  and  Kennedy  (1998);  and  Kuznetsov’s  equation  for  rock  blasting, 
e.g.,  Cunningham  (2005)  and  Ouchterlony  (2016). 
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